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ABSTRACT 

We present a study of the Lya forest at z < 0.4 from which we conclude that at least 20% of the 
total baryons in the universe are located in the highly-ionized gas traced by broad Lya absorbers. 
The cool photoionized low-z intergalactic medium (IGM) probed by narrow Lya absorbers contains 
about 30% of the baryons. We further find that the ratio of broad to narrow Lya absorbers is higher 
at z < 0.4 than at 1.5 < z < 3.6, implying that a larger fraction of the low redshift universe is hotter 
and/or more kinematically disturbed. We base these conclusions on an analysis of 7 QSOs observed 
with both FUSE and the fiST/STIS E140M ultraviolet echelle spectrograph. Our sample has 341 H I 
absorbers with a total unblocked redshift path of 2.064. The observed absorber population is complete 
for logiVni > 13.2, with a column density distribution /(iVm) oc N^f . For narrow (b < 40 kms -1 ) 
absorbers (3 = 1.76 ± 0.06. The distribution of the Doppler parameter b at low redshift implies two 
populations: narrow (b < 40 kms -1 ) and broad (b > 40 kms -1 ) Lya absorbers (referred to as NLAs 
and BLAs, respectively). Both the NLAs and some BLAs probe the cool (T ~ 10 4 K) photoionized 
IGM. The BLAs also probe the highly-ionized gas of the warm-hot IGM (T ~ 10 5 -10 6 K). The 
distribution of b has a more prominent high velocity tail at z < 0.4 than at 1.5 < z < 3.6, which 
results in median and mean b- values that are 15-30% higher at low z than at high z. The ratio of the 
number density of BLAs to NLAs at z < 0.4 is a factor of ~3 higher than at 1.5 < z < 3.6. 
Subject headings: cosmology: observations — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

Observations of Lya absorption lines in the spectra 
of QSOs provide a sensitive probe of the evolution and 
the distribution of the gas in the universe from high 
to low redshift. A forest of H I absorption lines oc- 
curs at different redshifts, z, along QSO sightlines with 
logA^Hi < 17. Understanding the evolution of the Lya 
forest with redshift is critical to understanding the evo- 
lution and formation of structures in the universe. At 
z > 1.5, observations of the Lya forest are obtained 
from ground-based telescopes at a spectral resolution of 
7-8 kms" 1 using 8- 10 m class telescopes (e.g ., iHu et all 
1995jlLu et aLlll996HKirkman fe Tvtlerl H997: Ki m et all 
19971 l2002aL but at z < 1.5 they require UV space- 



based instruments. Space-based UV astronomy has pro- 
duced remarkable results, including the discovery itself 
of the Lya forest at low redshift (B ahcall et al.l 119911 : 
iMorris eta l. 1991). However, most of the UV studies 
of the IGM at low z have lacked the spectral resolution 
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Ultraviolet Spectroscopic Explorer. FUSE is operated for NASA by 
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and wavelength covera ge of the hig her redshift studies 
(IWeymann et al.lll998t llmpev et al1ll999t iPenton et all 
I2000ll2004f) . requiring assumptions for the Doppler pa- 
rameter to derive the column density. The situation at 
z < 0.5 has dramatically improved in the last few years 
with high quality observations of several low redshift 
QSOs obtained with the Hubble Space Telescope (HST) 
and its Space Telescope Imaging Spectrograph (STIS). 
In its E140M echelle mode, STIS provides a spectral res- 
olution of ^7 kms -1 , comparable to the resolution of 
the high redshift observations of the Lya forest. The 
high spectral resolution has allowed the derivation of ac- 
curate Doppler parameters, b, and column densities, N, 
using techniques similar to those used at high redshift. 
The Doppler parameter is important because it is related 
to the temperature of the gas via b 2 = b 2 h + b 2 t , where 
6 nt is the non-thermal broadening of the absorption line 
and 6 t h = \/2kT/m = 0.129\/T is the thermal broad- 
ening of the H I absorption line; therefore the measured 
b directly provides an upper limit to the temperature of 
the observed gas. 

In parallel, hydrodynamic cosmological simulations of 
the local un i verse h ave quickly evolved in r ecent years 
dDave et al.l 119991: ICen fc Ostrikerl Il999t iDave et all 
12001b iCen fc Fand 120061 : iCen fc Ostriker! I2006T I . These 
models predict that at low redshift roughly 30-50% of 
the baryons are in a hot (10 5 -10 7 K) and highly ion- 
ized intergalactic medium (IGM) known as the warm-hot 
intergalactic medium (WHIM), 30-40% are in a cooler 
medium (< 10 4 K) photoionized by the UV background, 
and the remaining baryo ns are in galaxies. Observation- 
ally, IPenton et all (|2004f) found with moderate spectral 
resolution (FWHM~20 kms -1 ) UV observations that 
the cool phase of the IGM may contain 29% of the baryon 
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budget. The observational detection of the WHIM came 
first from an intensive search of collisionally-ionized O VI 
systems and other highly ionized species such as Ne VIII 
using space-based observatories such as the Far Ultravio- 
let Spectroscopic Explorer (FUSE) and HST/STIS (e g , 
Tripp. Savage, fc Jenkins! mM iDanforth fc STiul [20051 : 
Savage et alj l2005f T The baryon content of the O VI 



absorbers suggests that ilb(0 VI) > 0.0022/^70 or at 
least 5% of the total baryon budget (ISembach et all2004t 
iTripp et alJ 120061 : IDanforth fc Sh ull 2005), but this esti- 
mate relies critically on the assumed oxygen abundance 
and the number of collisionally io nized O VI systems ver- 
sus photoionized O VI systems (jProchaska et al.l 12004 
I Lehner et a l. 2006). The O VI systems probe the lower 
end of the WHIM temperature range (T < 10 6 K). 
Higher temperatures can be traced with more highly ion- 
ized oxygen ions (O VII, O VIII) that can be observed 
in p rinciple with X-ray observatories suc h as Chan- 
dra ([Nicastro et al.l 120051 : iFang et all I2006D and XMM- 
Newton, but the I GM detections of O VII and O VIII are 
still c ontroversial (Kaa stra et al.l [20061 iRasmussen et al.1 
[2001 for z > 0. 

The WHIM can also be detected through broad H I 
absorption lines. The high temperature of the WHIM 
will broaden the Lya absorption line resulting in a 
large Doppler parameter (b > 40 kms -1 ). A very 
small fraction of H I (typically < 10~ 5 ) is expected to 
be found in the highly ionized plasma of the WHIM, 
so that the broader the H I absorption line is, the 
shallower it should be. Recent observations, both 
with moderate spectral resolution and with the higher- 
resolution STIS E140M echelle mode, have revealed the 
presence of broad H I absorption lines that could be 
modeled by smooth and broad Gaussian components 
(iTripp et al.l200HlBowen et al.ll2002l:lRichter et all 2004: 
ISembach et al.l l2004t iLehner et al.ll2006h . H I absorbers 
with b < 40 kms" imply that the temperature of these 
absorbers must be T < 10 5 K. Since the nominal lower 
temperature of the WHIM is T ~ 10 5 K, it is a priori 
natural to consider two different physical populations of 
H I absorbers: the narrow Lya absorbers (NLAs) as H I 
systems having b < 40 kms~ and the broad Lya ab- 
sorbers (BLAs) as H I systems having b > 40 kms -1 . 
We will see that the separation between the NLAs and 
BLAs is not so evident observationally, with an impor- 
tant overlap between the two populations, particularly 
in the range b = 40-50 kms' 1 . While it is clear that 
the NLAs are tracing mostly the cool photoionized IGM, 
the situation is less clear for the BLAs. There is a 
fuzziness in the separation of the BLAs from the NLAs 
because non-thermal broadening can be important and 
unresolved components can hide the true structure of 
the H I absorbers. Recent simulations show that non- 
thermal broadening could be particularly important for 
the H I absorbers with 40 < b < 60 km s _1 ([Richter et al.l 
2006a). Therefore BLAs may probe photoionized gas 
(that can be either cool (T < 10 4 K) or hot (T > 10 5 K) 
but with a very low density (lognn < — 5.3)) and colli- 
sionally ionized, hot (T > 10 5 K) gas. In this paper, we 
define a BLA as an absorber that can be fitted with a 
single Gaussian component with b > 40 kms -1 (a NLA 
has b < 40 kms~ and can have multiple components). 
BLAs probe the IGM that is either hotter (6 t h 3* &nt) 



or is more kinematically disturbed (& n t 3> &th) than the 
IGM probed by the NLAs. Studying in detail the NLA 
and BLA populations is also important for estimating 
the baryon density because an accurate inventory of the 
baryon distribution must separate the fraction of the 
baryons that are located in the cool photoionized IGM 
versus those that are in the substantially hotter shock- 
heated WHIM phase. The narrow and broad Lya lines 
provide a means to discriminate between the cool and 
shock-heated gas clouds. Therefore to make a reliable 
assessment of the baryonic content of the Lya forest at 
low z, it is necessary to investigate the frequency and 
properties of the NLAs and BLAs. The current estimate 
of the baryon budget residing in the BLAs over a redshift 
path Az = 0.928 yields f} b (BL A) > 0.0027/t 7 " 11 1 or at least 
6% of the total baryon budget ([Richter et al.ll2006bh . as- 
suming that the observed broadening is mostly thermal 
and collisional ionization equilibrium applies. 

In this paper, we will address these issues using a sam- 
ple of 7 QSO sightlines that have been observed with 
both STIS E140M and FUSE and for which the data 
have been fully analyzed using similar techniques and 
are in press or to be submitted soon. While the spec- 
tral resolution of FUSE is only ~20 kms -1 compared to 
-7 kms -1 for STIS E140M, FUSE gives access to sev- 
eral Lyman series and metal lines making the line iden- 
tification more reliable and providing an unprecedented 
insight of the physical conditions and metallicity. The 
fundamental parameters of the Lya forest lines (redshift 
z, column density N, and Doppler parameter b) were 
accurately determined using profile fitting of all the ob- 
served Lyman series lines. Because b could be derived, 
this is the first time with a large sample that the prop- 
erties of the Lya forest in the low redshift universe can 
be studied as a function of b. The main aims of this 
work are 1) to study the distribution and evolution of 
the Doppler parameter of the Lya forest, and 2) to de- 
termine the baryon density of the Lya forest in the NLA 
and BLA populations. The organization of this paper is 
as follows: describes the sample and its completeness; 
in Sj3]we study the distribution of the Doppler parameter 
b and the Lya density number as a function of b; in S|4]we 
study the evolution of b with redshift by comparing the 
low-z (z < 0.4) Lya forest with the mid-z (0.5 < z < 1.5) 
and high-z (1.5 < z < 3.6) Lya forest; in Sj5]we estimate 
the column density distribution; and finally, in $6] we es- 
timate the baryon content of the cool photoionized IGM 
probed by the NLAs and the photoionized IGM and the 
WHIM probed by the BLAs, and discuss the uncertain- 
ties in estimating the baryon budget. We summarize our 
results in 

2. THE LOW Z SAMPLE 

2.1. Description of the Sample and Completeness 

Our low-z sample consists of 7 QSO sightlines that 
were observed with HST/STIS E140M and FUSE: 
H 1821+643 (Sembach et al. 2007, in preparation; 
see also Tripp, Savage, & Jenkins 2000 and Oegerle 
et al 2000 for the m etal-line systems ), HE 0226- 
Tnjj jLehner et al.l l20M ). HS 0624+6907 (|Aracil et all 
I2006aj) . PG 0953+415 (TM. Tripp 2006, private com- 
munication; see also Savage et al. 2002 for th e 
metal-line systems), PG 1259+593 (|Richter et all 120041 . 
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TABLE 1 
Line of Sight Properties 



Sightline 


S/N 


z 


Az 


AX 


Refs. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


H 1821+643 


15-20 


0.297 


0.238 


0.266 


(1) 


HE 0226-4110 


5-11 


0.495 


0.401 


0.481 


(2) 


HS 0624+6907 


8-12 


0.370 


0.329 


0.383 


(3) 


PG 0953+415 


7-11 


0.239 


0.202 


0.222 


('!) 


PG 1116+215 


10-15 


0.176 


0.126 


0.134 


(5) 


PG 1259+593 


9-17 


0.478 


0.355 


0.418 


(«) 


PKS 0405-123 


5-10 


0.574 


0.413 


0.498 


(7) 



Note. — Column 2: Range of signal-to-noisc (S/N) per 
resolution element in STIS E140M mode. Column 3: QSO 
redshift. Column 4: Unblocked rcdshift path for Lya. Col- 
umn 5: Absorption distance (see ^2j. Column 6: Refer- 
ences: (1) Scmbach ct al. (2007, in prep.), (2) Lehner et 
al. (2006), (3) Aracil ct al. (2006), (4) T.M. Tripp (2006, 
priv. comm.), (5) Scmbach ct al. (2004), (6) Richter et 
al. (2004), (7) Williger et al. (2006) and appendix of this 
paper. 



PG 1116+215 (|Sembach et al.ll200l . and PKS 0405-123 
(Williger et al. 2006; see also Prochaska et al. 2004 for 
the metal-line systems). Note that for HS 0624+6907, we 
us ed the result s summ arized in the erratum produced 
by lAracil et al.l ((2006b) . The data handling and anal- 
ysis are described in detail in the above papers. For 
PKS 0405-123, we adopt the new measurements and a 
new line list that we describe in the Appendix. The moti- 
vation to revisit Williger et al. 's analysis was first driven 
by differentiating a real detection from a noise feature 
for the BLAs since these authors noted that several of 
their BLAs could be just noise. This re-analysis also pro- 
vides an overall coherent data sample that was analyzed 
following the same methodology. Signal-to-noise (S/N) 
where Lya can be observed, redshift (z), unblocked red- 
shift path (Az), and absorption distance (AX) are sum- 
marized in Table [1] for each sightline. The absorption 
distance was computed assuming a Friedman cosmology, 
AX = 0.5[(1 + Az) 2 - 1] with q = 0. For the lines 
of sight to QSOs with zqso < 0.42, we have excluded 
absorption systems within 5000 kms -1 of the QSO red- 
shift. Some observations of low-redshift QSOs have pro- 
vided evidence that "intrinsic" absorption lines can arise 
in clouds that are spatially close to the QSO and yet, due 
to the cloud kinematics, are substantially offset in red- 
shift from the QSO (e.g., Yuan et al. 2002; Ganguly et 
al. 2003). To avoid contaminating the sample with these 
intrinsic absorbers, we do not use systems detected with 
5000 km s" 1 of the QSO. The Voigt profile fitting method 
was employed for each line of sight to measure the column 
densities (Nm), Doppler parameters (b), and redshifts of 
the absorbers and we adopt those results for our anal- 
ysis. Table H lists these parameters. For PG 1259+593 
we did not include in our sam ple the systems ma rked 
uncertain (UC) in Table 5 of iRichter etail (|2004l ). A 
colon in Table [5] indicates that there is uncertainty in 
the determination of the physical parameters, which are 
not accounted for in the formal errors produced by the 
profile fitting. These systems are not taken into account 
when we use an error cutoff (see below). If we had, it 
would not have changed the results presented in this pa- 
per in a statistically significant manner. 

The total sample consists of 341 H I systems, with a 
total unblocked redshift path EAz = 2.064 and a total 



absorption distance EAAT = 2.404. There are 201 sys- 
tems at z < 0.2 and 131 systems at 0.2 < z < 0.4. The 
remaining 9 systems lie at 0.4 < z < 0.44. Therefore, 
our sample mostly probes the universe at z < 0.4. More 
absorbers are found at z < 0.2 because several lines of 
sight do not extend out to z = 0.3-0.4 and to a lesser 
extent because the S/N at A > 1650 A decreases rapidly 
(see below). 

TABLE 2 

H I MEASUREMENTS 



z 


log A/hi 


b 




(dex) 


(kms" 1 ) 




H1821- 


h643 (Scmbach et al. 


2006, in prep.) 


0.02438 


14.28+°'°! 


26.9 + 1.7 


0.02642 


13.26 


48.6 + 6.0 


0.06718 


13-72 +g : gj 


20.2 + 0.6 


0.07166 


13.87 +° : <g 


34.0 + 1.0 


0.08911 


13.01 +^ 6 7 


23.3 + 2.8 


0.11133 


12.95 


88.0 + 14.0 


0.11166 


12-99 


28.6 + 4.4 


0.11961 


13.15 


36.8 + 3.4 


0.12055 


12.64 t°:ll 


22.5 ± 5.9 


0.12112 


13.93 + 0.37 


26.0+i 3 ° 


0.12125 


14.04 ± 0.36 


40.01&0 


0.12147 


13.78 ± 0.17 


85.0 ^26.0 


0.12221 


13.19 +°f 6 


41.7 + 4.4 


0.14754 


13-51 +0 f 3 


44.6 + 2.5 


0.14776 


13.30 +g-g* 


19.2 ± 1.6 


0.15731 


13.H ±g:Sg 


22.9 + 2.1 


0.16127 


12.74 ±g-i° 


21.5 + 4.0 


0.16352 


in i 7 +0.06 
ro.ri _q 07 


52.2 ± 5.9 


0.16966 


13-93 ±°-jg 


35.7 + 2.0 


0.17001 


13 65 + 02 
lO.DO _ .02 




0.17051 


1 O AT* +0.04 

13.40 ±o.o4 


21.1 ± 1.4 


0.17926 


12.87 ±°;» 


18.4 + 3.4 


0.18047 


i3-i4lS:£ 


50.7 + 6.6 


0.19662 


12-98 to 1 ;^ 


21.1 ± 2.7 


0.19904 


12 74 +0 ' 12 
lz - ' 4 -0.17 


22.1 ± 5.4 


0.20957 


,o ln +0.03 
Id.lU _ .04 


16.2 + 0.9 


0.21161 


13-16 ±°;° 6 6 


22.3 + 2.3 


0.21668 


12.88 


19.7 + 3.0 


0.21326 


14 41 +0 04 
14.41 -0.04 


43.0 + 2.4 


0.22497 


It KQ+0.05 
10. OO _ .05 


25.0 + 7.0 


0.22616 


la C] +0.04 


54.7 + 3.8 


0.22786 


13-26 ±°;°3 


35.3 + 2.0 


0.23869 


12-86 t°j° 


20.2 + 3.8 


0.24142 


io i 9 +0.03 
lO.lz _ .04 


23.4 + 1.4 


0.24531 


13-06 ±°;°g 


34.2 ± 5.1 


0.25689 


12.80 ±°;il 


22.6 + 5.1 


0.25814 


13-38 toTs 


60.3 + 8.8 


0.25816 


12-98 tgis 


14.5 ± 3.4 


0.26152 


,o 7n +0.02 
Id. fU _ .03 


37.7 + 1.5 


0.26659 


13 64 + 03 
lo.O4_ 03 


44.5 ± 2.1 


HE 0226-4110 (Lehner et al. 2006) 


0.01746 


13.22 ± 0.06 


17.9 + 4.3 


0.02679 


13.22 + 0.08 


41.6 + 11.0 


0.04121 


12.82 ± 0.14 


23.6 + 18.6 


0.04535 


12.71 ± 0.13 


16.2 ± 12.8 


0.04609 


13.66 + 0.03 


25.0 + 2.1 


0.06015 


13.19 ± 0.06 


35.5 ± 7.6 


0.06083 


14.65 ± 0.02 


44.5 ± 1.0 
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TABLE 2 — Continued 





z 


log TVhi 
(dex) 


b 

(kms -1 ) 





.07023 


13.81 ±0.11 


26.0 ± 12.4 





.08375 


13.67 ±0.05 


29.6 ±4.5 





.08901 


13.33 ±0.05 


23.8 ± 3.7 





.08938 


12.59 ±0.22 


8.2 : 





.08950 


12.72 ±0.20 


22.0 : 





.09059 


13.71 ±0.03 


28.3 ±2.0 





.09220 


12.94 ±0.11 


40.2 ± 18.0 





.10668 


13.09 ±0.08 


32.7 ±9.2 





.11514 


12.90 ±0.09 


10.4 ±4.1 





.11680 


13.27 ±0.05 


23.7 ±3.9 





.11733 


12.64 ±0.15 


15.0 : 





.12589 


13.01 ±0.09 


29.2 ± 10.1 





.13832 


13.19 ±0.06 


25.9 ±5.3 





.15175 


13.42 ± 0.05 


48.6 ± 6.7 





.15549 


13.13 ±0.08 


34.7 ±9.8 





.16237 


13.04 ±0.08 


29.7 ±8.6 





.16339 


14.36 ±0.04 


46.3 ± 1.9 





.16971 


13.35 ±0.05 


25.3 ±3.9 





.18619 


13.26 ±0.08 


53.9 ± 16.2 





.18811 


13.47 ± 0.05 


22.4 ±3.3 





.18891 


13.34 ±0.07 


22.2 ±4.0 





.19374 


13.20 ±0.06 


28.7 ±6.0 





.19453 


12.89 ±0.12 


26.1 ± 14.0 





.19860 


14.18 ±0.04 


37.0 ±2.0 





.20055 


13.38 ±0.05 


38.9 ± 6.4 





.20698 


13.31 ±0.34 


97.0 : 





.20700 


15.06 ±0.04 


17.4 ± 1.4 





.20703 


14.89 ± 0.05 


35.9 ± 1.1 





.22005 


14.40 ± 0.04 


27.7 ± 1.1 





.22099 


12.99 ±0.12 


34.1 ± 18.1 





.23009 


13.69 ±0.04 


67.9 ±7.5 





.23964 


13.13 ±0.08 


28.8 ± 8.8 





.24514 


14.20 ±0.03 


34.5 ± 1.6 





.25099 


13.17 ±0.08 


37.9 ± 11.5 





.27147 


13.85 ±0.07 


25.7 ±4.2 





.27164 


13.33 ±0.28 


26.2 : 





.27175 


12.88 ±0.38 


11.2 : 





.27956 


13.22 ±0.14 


36.3 ± 24.6 





.28041 


13.03 ±0.11 


13.9 ± 6.7 





.29134 


13.53 ±0.07 


27.0 ±6.2 





.29213 


13.19 ±0.12 


33.4 ± 17.8 





.30930 


14.26 ±0.03 


43.8 ± 2.3 





.34034 


13.68 ± 0.06 


33.4 ± 4.9 





.35523 


13.60 ±0.07 


27.1 ± 6.8 





.37281 


13.16 ±0.12 


25.9 ± 13.3 





.38420 


13.91 ± 0.04 


62.0 ±7.1 





.38636 


13.36 ±0.09 


38.1 ± 12.7 





.39641 


13.59 ±0.10 


62.8 ±22.7 





.39890 


13.50 ±0.16 


151.7 : 





.40034 


13.39 ±0.11 


60.7 ±26.1 





.40274 


14.13 ±0.04 


45.7 ±4.2 


HS 0624+6907 (Aracil et al. 2006b 1 





.01755 


12.96 ±0.05 


29.0 ±4.3 





.03065 


13.36 ±0.03 


22.0 ± 1.7 





.04116 


13.33 ±0.03 


41.0 ±3.0 





.05394 


13.26 ±0.04 


24.0 ±2.3 





.05437 


13.09 ±0.11 


60.0 ±19.2 : 





.05483 


14.50 : 


35.0 : 





.05515 


13.68 ±0.17 


84.0 ± 30.7 : 





.06188 


13.77 ±0.03 


21.0 ±1.4 





.06201 


12.63 ±0.17 


8.0 ±4.7 





.06215 


12.41 ±0.22 


10.0 ±7.9 





.06234 


13.45 ±0.05 


30.0 ±4.0 





.06265 


13.31 ±0.14 


35.0 ± 12.3 





.06276 


12.95 ±0.28 


8.0 ±3.7 





.06285 


13.42 ±0.14 


20.0 ±7.0 





.06304 


13.33 ±0.13 


27.0 ±8.8 





.06346 


14.46 ± 0.30 


48.0 ±8.4 : 





.06348 


15.27 ±0.13 


24.0 ±5.5 





.06362 


14.29 ±0.38 


10.0 ±5.6 





.06475 


13.87 ±0.04 


33.0 ±3.0 





.06502 


13.97 ±0.04 


31.0 ±2.7 





.07573 


14.18 ±0.03 


24.0 ± 0.8 





.09023 


13.29 ±0.08 


76.0 ± 13.7 
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z 


log Nui 


b 






(dex) 


(kms -1 ) 





.13076 


13.34 ±0.04 


34.0 ±3.6 





.13597 


13.33 ±0.10 


57.0 ± 10.7 





.16054 


13.08 ±0.21 


34.0 ± 10.3 





.16074 


13.66 ±0.05 


30.0 ±2.4 





.19975 


13.24 ± 0.05 


17.0 ±2.0 





.19995 


13.17 ±0.06 


26.0 ±4.6 





.20483 


13.72 ± 0.02 


24.0 ± 1.0 





.20533 


14.12 ±0.03 


25.0 ±0.8 





.20754 


13.48 ±0.02 


27.0 ± 1.5 





.21323 


13.22 ± 0.05 


45.0 ±5.6 





.21990 


13.39 ±0.05 


60.0 ±8.6 





.22329 


13.86 ±0.02 


25.0 ±0.9 





.23231 


13.33 ±0.08 


44.0 ±7.7: 


0. 


.23255 


12.86 ± 0.21 


24.0 ±7.3 


0. 


.24060 


13.33 ±0.04 


20.0 ±2.0 





.25225 


12.96 ±0.06 


24.0 ±4.2 





.26856 


13.03 ±0.05 


51.0 ±7.2 





.27224 


12.80 ±0.06 


12.0 ±2.2 





.27977 


13.50 ±0.06 


34.0 ±4.9 





.28017 


14.32 ± 0.02 


43.0 ± 1.9 





.29531 


13.80 ± 0.02 


42.0 ±2.0 : 





.29661 


13.54 ±0.02 


52.0 ±2.9 





.30899 


13.49 ± 0.03 


28.0 ± 1.8 





.30991 


13.61 ±0.10 


66.0 ±12.3 : 





.31045 


13.43 ± 0.33 


62.0 ±40.3 : 


0. 


.31088 


13.13 ±0.43 


51.0 ± 27.7 





.31280 


13.65 ± 0.10 


54.0 ±9.3 





.31303 


13.09 ± 0.24 


17.0 ±6.8 





.31326 


13.62 ± 0.10 


55.0 ± 10.9 





.31790 


13.37 ±0.04 


34.0 ±3.6 





.32089 


13.97 ±0.02 


31.0 ± 1.2 





.32724 


13.73 ±0.32 


69.0 ±15.6 : 





.32772 


13.61 ± 0.43 


115.0 ±62.1 : 





.33267 


13.55 ± 0.04 


38.0 ± 3.4 


0. 


.33976 


14.45 ±0.03 


42.0 ± 1.3 


0. 


.34682 


13.59 ±0.02 


39.0 ± 1.9 





.34865 


12.78 ± 0.06 


18.0 ±3.0 



PG0953±415 (T.M. Tripp, 2006, priv. comm.) 






.01558 


13 


.19 


± 


0. 


08 


35.0 ±iy 

19.0 ±1™ 





.01587 


12 


.86 


± 


0. 


.21 





.01606 


13 


.51 


± 


0. 


05 


32.0 ±5.0 





.01655 


13 


.51 


± 


0. 


.04 


26.0 ±3.0 





.02336 


13 


.21 


± 


0. 


.08 


56.0 ±»;° 





.04416 


12 


.72 


± 


0. 


.09 


15.0 ±6.0 





.04469 


12 


.77 


± 


0. 


06 


12.0 ±3.0 





.04512 


13 


.42 


± 


0. 


02 


38.0 ±2.0 





.05876 


13 


.82 


± 


0. 


07 


25.0 ±3.0 





.05879 


13 


.41 


± 


0. 


16 


63.0 ±»;g 





.06808 


14 


.47 


± 


0. 


.03 


21.0 ±2.0 





.09228 


13 


.08 


± 


0. 


.07 


27.0 ±6.0 





.09315 


13 


.66 


± 


0. 


.03 


39.0 ±3.0 





.10940 


13 


.68 


± 


0. 


.03 


23.0 ±2.0 





.11558 


13 


.47 


± 





.03 


24.0 ±2.0 





.11826 


13 


.68 


± 


0. 


02 


30.0 ±2.0 





.11871 


12 


.81 


± 


0. 


10 


16.0 ±6.0 





.12558 


12 


.77 


± 


0. 


.09 


s n+ 4 - 

S - U -3.0 





.12784 


12 


.83 


± 


0. 


.35 


44 n +76.0 
44. U _ 28 .o 





.12804 


13 


.27 


± 


0. 


12 


21.0 ±5.0 





.14178 


12 


.68 


± 





10 


10.0 1*;° 





.14233 


13 


.58 


± 


0. 


.03 


28.0 ±2.0 





.14263 


13 


.45 


± 





.04 


31.0 ±4.0 





.14294 


12 


.63 


± 





27 


17 0+ 23 

L '- u -10.0 





.14310 


13 


.05 


± 


0. 


.11 


22 +9.0 





.14333 


12 


.77 


± 





.14 


i5.o tn 





.17985 


13 


.27 


± 


0. 


07 


48.0 ±ft° 





.19072 


13 


.01 


± 





08 


26.0 ±6.0 





.19126 


13 


.08 


± 


0. 


.55 


48.0 ±l|;g 





.19147 


13 


.33 


± 


0. 


.34 


30.0 ±l 2 ° 





.19210 


13 


.11 


± 


0. 


07 


28.0 ±6.0 





.19241 


12 


.91 


± 


0. 


.12 


35.0lif;g 
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TABLE 2 — Continued 



z 


log 7V H i 


b 




(dex) 


(kms^ 1 ) 


0.19361 


13.94 + 0.02 


40.0 + 2.0 


0.20007 


13.24 + 0.09 


66 0+ 17 
DDU -14.0 


0.20104 


13.16 + 0.16 


71 +41.0 
'1-U -26.0 


0.20136 


12.91 + 0.19 


19.0+ 7 ° ° 


0.20895 


12.94 + 0.10 


28.0 jf ° 


0.21514 


13.30 + 0.04 


27.0 + 3.0 


0.22526 


12.74 + 0.23 


2 0+ 2 
z - u -1.0 


0.22527 


13.19 + 0.05 


34.0 + 5.0 


PG 1116+215 fSembach et al. 2004) 


0.00493 


13.36 

j-o.ou —0.06 


34.2 ± 3.6 


0.01635 


13.39 + 0.06 


48.5 ± 5.1 


0.02827 


13.80 + 0.02 


31.4 + 1.1 


0.03223 


10 oq+0.06 
lO.OO _o.05 


31.6 ± 2.9 


0.04125 


13 25 +0 11 


105.0 + 18.0 


0.04996 


10 7o +0.10 
1Z.(Z _q 8 


16.5 ± 3.2 


0.05895 


13.56 + 0.05 


25.0 + 5.0 


0.05928 


-0.13 


10.0 : 


0.06072 


i, Qo+0.06 


55.4 + 5.8 


0.06244 


13-18 _ 06 


77.3 ± 9.0 


0.07188 


12 7q +0.08 
iz. in _ .07 


9.6 + 2.0 


0.08096 


13 45 +0.03 


24.9 + 1.0 


0.08587 


12 90+ 19 


52.0 + 14.0 


0.08632 


1 2 66 +°- 32 


36.0 + 15.0 


0.09279 


13 no +0-09 


133.0 + 17.0 


0.10003 


19 7 , +0.07 

1Z. IO _Q Qg 


23.2 ± 2.2 


0.11895 


13 44+0-04 


31.5 ± 1.8 


0.13151 


13.41 + 0.03 


28.8 ± 1.3 


0.13370 


io 97+0.08 
10 ' z ' -0.07 


83.6 ± 10.4 


0.13847 


16 20 + 05 
10. ZU _ 04 


22.4 + 0.3 


PG 1259+593 fRichter et al. 20041 


0.00229 


13.57 + 0.10 


42.1 + 4.4 


0.00760 


14.05 + 0.05 


34.6 ± 2.0 


0.01502 


13.21 + 0.06 


22.6 + 4.4 


0.02217 


13.67 + 0.04 


30.2 ± 2.2 


0.03924 


12.94 + 0.05 


15.3 ± 2.8 


0.04606 


15.58 + 0.21 


47.6 + 12.4 


0.05112 


13.62 + 0.07 


34.4 + 2.5 


0.05257 


12.75 + 0.06 


20.7 + 3.0 


0.05376 


13.44 + 0.04 


30.5 ± 1.9 


0.06644 


13.65 + 0.05 


28.7 + 3.3 


0.08041 


12.97 + 0.10 


42.0 + 4.5 


0.08933 


14.04 + 0.03 


28.9 ± 1.7 


0.09591 


12.97 + 0.03 


21.5 ± 2.3 


0.12188 


13.03 + 0.07 


26.9 + 4.2 


0.12387 


13.47 + 0.06 


28.2 ± 3.0 


0.14852 


13.91 + 0.06 


42.1 ± 2.4 


0.15029 


13.25 + 0.11 


25.7 + 4.3 


0.15058 


13.45 + 0.13 


32.0 + 5.1 


0.15136 


13.32 + 0.09 


65.3 ± 5.5 


0.15435 


13.22 + 0.04 


25.2 ± 1.9 


0.17891 


13.29 + 0.10 


98.5 ± 9.1 


0.18650 


13.02 + 0.03 


19.6 ± 1.4 


0.19620 


13.65 + 0.05 


32.7 + 3.2 


0.19775 


13.33 + 0.06 


23.9 ± 2.8 


0.21949 


15.08 + 0.08 


32.3 ± 1.4 


0.22313 


13.92 + 0.04 


34.8 ± 1.1 


0.22471 


13.59 + 0.06 


28.9 ± 1.7 


0.22861 


13.47 + 0.05 


40.3 ± 2.9 


0.23280 


13.50 + 0.07 


37.4 + 3.2 


0.23951 


13.40 + 0.03 


16.3 ± 1.3 


0.24126 


13.41+0.09 


89.1 ± 6.9 


0.25642 


13.74 + 0.03 


25.0 + 0.9 


0.25971 


13.84 + 0.12 


40.5 + 4.9 


0.28335 


13.59 + 0.10 


37.0 + 5.2 


0.29236 


14.65 ± 0.09 


24.3 ± 2.8 


0.29847 


13.09 + 0.10 


33.3 ± 2.9 


0.30164 


13.26 + 0.14 


31.7 + 4.7 



TABLE 2 — Continued 





z 


log Nm 


b 






(dex) 


(kms -1 ) 





.30434 


13.76 ± 0.14 


64.5 + 9.6 





.31070 


13.40 ± 0.07 


22.8 + 2.9 





.31978 


13.98 ± 0.06 


74.4 + 8.7 





.32478 


13.24 + 0.15 


46.1 ± 10.2 





.33269 


13.88 + 0.08 


25.9 + 3.2 





.34477 


14.02 ± 0.08 


34.3 ± 4.4 





.34914 


13.36 ± 0.09 


31.3 + 4.8 


Q 


.35375 


13.41 ± 0.06 


16.4 + 1.8 





.37660 


13.45 ± 0.13 


36.4 + 6.2 





.38833 


13.02 + 0.19 


14.1 + 3.8 





.41081 


13.57 + 0.08 


32.8 + 5.0 





.41786 


13.25 ± 0.08 


50.7 + 3.9 





.43148 


14.10 + 0.06 


20.9 + 4.3 





.43569 


14.22 ± 0.10 


44.0 + 3.9 




PKS 0405-123 (this paper, 


see Appendix) 



The sample is not homogeneous with respect to the 
achieved S/N in the STIS E140M observations toward 
the 7 sightlines considered. The detection limit depends 
on the S/N and the breadth over which the spectrum is 
integrated. The matter is complicated by the S/N not 
being constant over the full wavelength range from ^1216 
A to ~ 1730 A available with STIS E140M where Lya 
can be observed; in particular, it deteriorates rapidly at 
A > 1650 A. Only lines of sight with z > 0.32 can reach 
the last 100 A of the STIS E140M wavelength coverage. 
Note that 297 systems with z < 0.32 are observed at 
wavelengths A < 1650 A. For example, in one of the 
lowest S/N spectra in our sample (HE 0226-4110), we 
estimate that at ~1700 A, a 3<t limit is ~75 mA for the 
profile integrated over Sv = [—50,50] kms -1 and ~100 
mA over 5v = [—90,90] kms -1 . We estimate that our 
sample is complete for log iV(rI I) > 13.20 (corresponding 
to a rest frame equivalent width W1215 — 88 mA) at 
a 3cr level for b < 80 kms -1 . For the high S/N lines 
of sight (H 1821+643, PG 1259+593, and PG 1116+215), 
this limit is quite conservative. For systems with b > 80 
kms -1 and log 7V(H I) > 13.20, our sample is incomplete, 
especially at z > 0.32 and for the lowest S/N spectra. 

Since our sample is not homogeneous with respect to 
the achieved S/N, it is useful to have a sample in which 
the cloud parameters are relatively well determined so 
that scatter due to noise is reduced. We therefore con- 
sider systems with errors on b and ./Vhi that are less than 
40%. The sample has 270 H I systems with errors on b 
and JVhi that are less than 40%. In this case, there are 
155 systems at z < 0.2, 107 systems at 0.2 < z < 0.4, 
and 8 systems at 0.4 < z < 0.44. At the completeness 
level log7V(H I) > 13.20 and with errors on b and Nm 
less than 40%, there is a total of 202 H I systems with 
109 systems at z < 0.2, 85 systems at 0.2 < z < 0.4, and 
8 systems at 0.4 < z < 0.44. 

2.2. Overview of the Distributions of Nm> b, z 

In Fig. [TJ we show the distribution of the column den- 
sity for each line of sight and the total sample (last panel, 
where we show the completeness limit of the sample) for 
systems with N m < 10 14 cm -2 and a b /b,a N /N < 0.4. 
About 86% of the absorbers have Nm < 10 14 cm -2 and 
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1. — Distribution of the column density of H I for each sightline and the combined sample with a column density interval of 6 X 10 12 



cm -2 . Only data with N m < W 1 



sample marks the completeness limit of Nm = 1.6 X 10 13 

~94% of them have N m < 10 14 5 cm- 2 . The column 
density distribution peaks near JVhi ~ 2 x 10 13 cm -2 
and drops sharply at smaller column densities. Since the 
peak of the distribution corresponds to about our com- 
pleteness limit, the observed decrease of the number of 
systems at TVhi < 1.6 x 10 13 cm -2 can be understood 
from the reduction in sensitivity. Higher S/N spectra 
would be needed to further understand the distribution 
of the smaller column densities. 

In Fig. [2j we show the distribution of the Doppler pa- 
rameters for each line of sight and for the entire sample 
(last panel) with <j b /b 7 <j N /N < 0.4. For every sight- 
line, the ^-distribution peaks between about 20 and 30 
kins -1 . Yet, it is apparent that the distribution is not 
Gaussian around these values, and, in particular, there 
are many systems with b > 40-50 kms -1 that produce a 
tail in the b distribution. This is clearly observed in the 
histogram of the combined sample where b peaks around 
20-30 kms -1 with an excess of systems with b > 40 
kms -1 compared to the number of systems with small b. 
This preview already shows clearly that a mean 6-value 
does not provide an adequate description of the Lya for- 
est. 

In Fig. [31 we show the distribution of the redshift for 
each line of sight with a redshift bin of 0.0033 (corre- 
sponding to 1000 kms - intervals). There is no strik- 
ing difference between the different lines of sight, but 
there are clearly voids and clustering in the distribu- 
tion, consistent with the Lya forest tracing filaments of 
matter in the universe. An example of clustering of ab- 
sorbers is found near the system at z — 0.06352 toward 



and less than 40% errors in b and Nm are considered. The vertical dashed line on the complete 



HS 0624+6907 where at least 8 absorbers are found in a 
single bin an d 17 absorbers are on ly separated by about 
3000 kms" 1 . lAracil et all (|2006al ) attribute this cluster- 
ing to the absorption of intragroup gas, possibly from a 
filament viewed along its long axis. 

3. DISTRIBUTION OF THE DOPPLER PARAMETER 

Earlier low-redshift UV studies of the H I forest have 
low or moderate spectral resolutions and smaller wave- 
length coverage and did not allow access to several Ly- 
man series transitions. Therefore, the Doppler parameter 
b generally had to be assumed and a study of the evo- 
lution and distributi on of b was not possible (see for ex- 
ample the s tudies of lWevmann et al.lll998t iPenton et al.1 
[20001 l200l . With STIS E140M observations (spectral 
resolution of 6.5 kms -1 ), b can be derived from profile 
fitting analysis. Furthermore, combining STIS E140M 
and FUSE observations allows further constraints on b 
by using several Lyman series lines and reducing possi- 
ble misidentifications. Here, we review the frequency and 
properties of the narrow Lya absorbers (NLAs, b < 40 
kms -1 ) and the broad Lya absorbers (BLAs, b > 40 
kms -1 ). 

3.1. Distribution of b and Other Low z Studies 

The median b- value of 31 km s -1 for our combined sam- 
ple is larger than the medians f o und in the l ow redshift 
IGM st udies bv iDave fe Trfpd (|2001h and IShull et all 
(2000). IDave k. Tr ipp (200^used automated software 
to derive b and N, but their criteria did no t allow a search 
for broad components. IShull et"aTI ([2000 < ) went only after 
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60 80 100 120 140 20 40 60 80 100 120 140 

b (km s-) b (km s -1 ) 

Fig. 2. — Distribution of the Doppler parameter of H I for each sightline and the combined sample. The observations arc plotted in 6 
kms -1 bins. Only data with less than 40% errors in b and JVhi are considered. 
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Fig. 3. — Distribution of the redshift of H I for each sightline with a redshift interval of Az = 0.0033. The arrows show the rcdshifts of 
the QSOs. The wavelength coverage of the E140M mode of STIS limits Lya measurements to z < 0.42. 
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log N(H I): Q y 12A0 

■ > 13.00 

■ > 13.20 

■ > 13.64 



40 60 80 100 120 140 

Doppler parameter (km s~ L ) 

Fig. 4. — Distribution of the Doppler parameter of H I for a 
series of minimum logA r (H I) for z < 0.4 for the seven sightlines. 
Only data with less than 40% errors in b and Nm are considered. 

the Ly/3 absorption lines in the FUSE wavelength range 
to combine with known Lya absorption lines observed 
with GHRS and we re therefore less like ly to find the 
broad H I absorbers. iDave fc Tripd (|2001[) found median 
and mean b- values of 22 kms -1 and 25 kms -1 , respec- 
tively. We find that the distribution of b is not gaussian, 
making the mean and dispersion less useful quantities. If 
we restrict our sample to data with b < 40 kms -1 , the 
median and mean (with la dispersion) are 27 kms -1 and 
27 ± 7 kms -1 . If only data with log A(H I) > 13.2 are 
considered, we obtain 29 kms" 1 and 28 ± 6 kms -1 . If 
we set the cutoff at b < 50 kms -1 , the median and the 
mean increase by 2 kms -1 . The median, mean, and dis- 
persion of b in our sample (with b < 40 or 5 km s -1 ) 
compare well to those derived bv lShull et al.l (|2000f k 28 
kms - and 31 ± 7 kms - for the median and mean, re- 
spectively. 

In Fig. IU we show the distribution of b for samples with 
various Ahi cutoffs. For any JVhi cutoff, the maximum of 
the distribution always peaks near 25-30 kms -1 and in 
all cases there is clearly an asymmetry in the distribution 
with the presence of a tail in the distribution that devel- 
ops at b > 40-60 kms -1 . For the BLAs, the number of 
systems with H I column between 13.0 and 13.2 dex is 
the largest, although this could be in part an observa- 
tional bias since weaker column density absorbers with 
b > 40 kms -1 would require higher S/N to be detected. 
Although the effects of line blending can contaminate the 
measurements of the observe d tails, the various p apers 
describing the data (see also lRichter et al.ll2006bf ) show 
that a large fraction of these broad absorbers can be de- 
scribed with a single Gaussian within the S/N. For the 
stronger of these absorbers, several Lyman series lines 
were also used to derive the physical parameters. 

3.2. The b-Nni distribution 

The b-Nm distribution of the low redshift sample is 
shown in Fig. [5l The solid gray line shows the threshold 
detection of H I absorbers. This curve is the relation- 
ship between b and Ntn for a Gaussian line profile with 



10% central optical depth - absorption lines to the left 
of this line are not detectable. The lack of data with low 
b and fog Am < 12.7 is most likely because our sample 
is not complete at these low column densities. When 
absorbers with log Am > 13.2 and b > kms -1 are 
considered, the fr-Am plot reveals mostly a scatter dia- 
gram. Most of the absorbers are present in the column 
density range 13.2-14 dex. In Table [3j we list the me- 
dian, mean, and dispersion of b for the entire sample, 
NLAs, and BLAs. For the entire sample, as expected 
from Fig. [5j it is not clear if b increases or decreases as 
Am increases. A Spearman rank-order correlation test 
on the entire sample with log Am > 13.2 shows a very 
marginal negative correlation between b and Am with a 
rank-order correlation coefficient r = —0.08 and a sta- 
tistical significance t — 0.25. We note that when sys- 
tems with log Am ^ 13.2 are considered, it creates an 
apparent correlation (r = 0.21 and t — 4.7 x 10~ 4 ) be- 
tween Am and b, which can be understood in terms of 
measurement biases, since weak broad systems are more 
difficult to detect than weak narrow absorbers. There is 
no clear separation between the NLAs and BLAs in this 
figure, although we note that most BLAs are found at 
logA H i < 14.0 and no BLAs with b > 50 kms -1 have 
log Am > 14.0 (see below). This scatter and absence of 
clear separation between NLAs and BLAs are expected 
if the H I lines trace systems with different temperatures 
and turbulent velocities. 

For NLAs with b < 40 kms -1 and logA H i > 13.2, we 
show in Fig. O the median, mean, and dispersion of b 
(red curves and symbols) derived in six intervals of Am . 
These estimates are summarized in Table |3] This shows 
evidence of an increase of b with increasing Ahi for the 
NLAs, at least for the weak absorbers with log Am < 
14.1. The Spearman rank-order correlation test for the 
NLA sample shows a weak correlation between b and Am 
with r = 0.12 and t = 0.15. Since the sample is complete 
for the NLAs, the increase of b with increasing Am must 
be real. The large scatter is again expected if the H I 
lines are broadened as a result of different temperatures 
and turbulent velocities. 

In Fig. we present the 6-Ahi distribution only for 
the BLAs (left panel) , which shows that b appears to de- 
crease with increasing A/hi: For 13.1 < log Am < 13.5 
(this range is highlighted in Fig. [6] by the vertical dotted 
lines), b is distributed between about 40 and 130 kms - ; 
for 13.5 < log Am < 14.0, b is mostly distributed be- 
tween about 40 and 80 kms -1 ; and for log Am > 14.0, 
b is always lower than 50 kms" 1 . This trend is also 
confirmed in the last two columns of Table [3] (note that 
in the [13.8,14.1] interval there are only 6 systems, with 
3 of them having b > 60 kms - and the other 3 hav- 
ing b ~ 40 kms -1 ). The Spearman rank-order correla- 
tion for the BLA sample with log Am > 13.2 confirms a 
negative correlation between b and Am (r = —0.30 and 
t = 0.016). 

In the right panel of Fig. [6l we show the r e cent si mu- 
lation of BLAs undertaken bv iRichter et al.l (|2006al ). in 
which artificial spectra were generated from the hydrody- 
namical simulation. Their numerical model was part of 
an earlier investig ation of the O VI ab sorption arising in 
WHIM filaments (|Fang fc Brvanll200ff ). and they include 
collisional ionization and photoionization processes. The 




Fig. 5. — Line width versus H I column density for the Lya absorbers along the seven sightlines. Only data with less than 40% errors 
in b and iVfji are considered. The solid gray line shows b vs N(R I) for a Gaussian line shape with 10% central optical depth - absorption 
lines to the left of this line a re not detectable. The dotted gray line shows the predicted minimum b- value from hydrodynamic cosmological 
simulations at low redshift (Dave ct al. 199 3). The dot-dashed and long- dash ed lines shows the fits to lower cutoff in the b distribution 
of high redshift systems (Kirkman & Tytler 1997; Kim ct al. 2002b) (sec i|3,2l for more details). The red histogram and red symbols show 
the median b (triangles), the mean and dispersion of b (circles and vertical bars) in six bins of JVhi only for data with b < 40 kms -1 (see 
Table 0. 



simulated sample presented in Fig. [^corresponds to their 
high-quality sample that includes 321 BLAs with almost 
perfect Gaussian profiles (note that 58% of the sample 
has logiVni < 13.2). Since our observations are com- 
plete only to logiVni > 13.2, it is not surprising that 
we are missing absorbers below this limit. It is, how- 
ever, interesting to note that the simulation and obser- 
vations have a similar trend: (i) most of the broadest 
absorbers are found at logTVui < 13.5, and (ii) most of 
the strong absorbers (log Ahi > 14) have b < 50 kms -1 
(although we note that in the simulation a few systems 
have b up to 65 km s _ , and the simulation does not pro- 
duce BLAs with log-ZVui > 14.4). This is also in general 
ag reement with the simulation of the WHIM produced 
bv lDave et ail (|200lD where they show the WHIM frac- 
tion peaks for an overdensity p/p of ^5-30. If Eq. [2] 
(see below) applies for the BLAs, the H I column den- 
sity range 13.2 to 14.0 dex corresponds to p/p ~ 5-17, 
in gene ral agreement with the hydrodynamical simula- 
tions of iDave et al.l (|2001[ ) if the BLAs tra ce mostly the 
WHIM . We note that the simulation of iRichter et al.l 
(2006a) generally produces larger b than currently ob- 
served: for systems with logiVni > 13.2, the median, 
mean, and standard deviation are 59, 69, 31 kms -1 for 



the simulation, while they are 52, 57, 18 kms" 1 for the 
observations. The current S/N of the observed data lim- 
its the detection of the broader absorbers with b > 80 
kms -1 , as discussed in $2] We also note that the broader 
systems (6 > 80 kms -1 ) are likely to be more uncertain, 
especially since they are detected in the lowest column 
d ensity range ( l og iVm < 13.4). 

Sch ave et al.l (|1999f ) and others demonstrated that the 
temperature-density relation, which is well described by 
a power law T — Tq(p/ p) 7-1 (where p is the average den- 
sity, p/p is the overdensity of the IGM), implies a lower 
envelope to the b-Nm distribution. This lower envelope 
is not clearly observed in Fig. [5] We further explore this 
by considering the numerical cosmological simulations of 
the low-redshift Lya forest that predict that the temper- 
ature T4 = T/(10 4 K) is a power law of the overdensity 
p/p, with 

T 4 « 0.5 (J J , (1) 

for the coolest systems at any given density. The over- 
density is connected to the H I column density through 
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FlG. 6. — Left: Line width versus H I column density for the BLAs (fc > 40 kms -1 ) along the seven sightlincs. Only data with less 
th an 40% errors i n b and iVfji are considered. Right: Line width versus H I column density for the high-quality sample of BLAs simulated 
bv lRichter et all fl2006aD . The dashed vertical lines highlight the column density range [13.1,13.5] dex where the broadest absorbers with 
b > 80 kms -1 are detected in the observational sample. 



(pave et al.lll999f) . 
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If we combine these two equations, we have: 



HI 



0.42 



10" 



(2) 



(3) 



The pure thermal Doppler broadening for H is 6th = 
0.129v / T\ so for photoionized hydrogen absorbers we can 
write 



0th 



22 



HI 



10 



-0.122 



kms 



(4) 



This latter relation is shown in Fig [5] with the dotted 
line where we set z = 0.2, which is about the mean and 
median z in our sample. Since the p—T fit to photoion- 
ized absorbers is for the coolest systems, this relation 
should provide a lower envelope to the b-Nm distribu- 
tion. The lower envelope to the observed b—Nm distri- 
bution roughly agrees with the numerical simulations, at 
least as long as logAni < 14.2. But the low redshift 
sample does not provide yet as sharp a lower envelope to 
the b-Nm d istribution as high redsh i ft samples do (see 
for example iKirkman fc Tvtlerl [l997t iKim et al.ll2002aft 
because there are still too few systems at the complete- 
ness level. In Fig. we show with the dot-dashed line 
the relatio n 6 min = 20 + 41o g [iVHi/ (10 14 cm~ 2 )1 kms" 1 
found by IKirkman fc Tvtlerl ([1997) at z = 2.7 and 
with the long dashed line the relation log6 m i n = 1.3 + 
0.090 lo g[jV H T/(10 14 cm" 2 )1 kms" 1 found bv IKim et all 
(|2002bh at z = 2.1 (smoothed power-law fit). These 
power laws provide a good approximation to the observed 
lower envelope of b at high redshift. At low redshift at 
logA^Hi > 13.2, a few absorbers lie below these fits, es- 
pecially at logA^Hi > 15. For systems with logiVHi < 15 
it is not clear if the lower b cutoff evolves with redshift. 



Finally, following [Dave fc Tripd (|2001h . if we compare 
the median in the various intervals with b < 40 kms - 
(red curve in Fig. [5]) and 6 t h defined in Eq. [H we find 
that 6th ~ (0.6 — 0.7)6 o bs for a typical absorber with 
logA^Hi < 14.2. Therefore, the contribution from ther- 
mal broadening is substantial for the low redshift NLAs. 
However, if the some BLAs actually trace cool photoion- 
ized gas, the non-thermal broadening will be dominant 
in these absorbers. 

3.3. Lya Line Density 

This sample provides the first opportunity to investi- 
gate the relative number of systems as a function of the 
Doppler parameter in the \ow-z IGM. In TableHJ we sum- 
marize dM /dz for each sightline and for the combined 
sample where our sub-samples have either b < 40 kms - 
or b < 150 kms -1 . For both 6-samples, we choose three 
different column density ranges: (a) [13.2, 14.0] dex, (b) 
[13.2, 16.5] dex and (c) [13.64, 16.5] dex. The lower limit 
of samples (a) and (b) corresponds to our threshold of 
completeness. The largest observed column density in 
our sample is about 10 16 5 cm" 2 . Hence sample (b) cor- 
responds to the combined sample with W > 90 mA, 
while sample (a) only covers the weaker Lya lines which 
may evolve differently than the stronger lines since the 
weaker lines may arise from tenuous gas in the IGM while 
the stronger lines may mostly trace the gas in the out- 
skirts of galaxies. The threshold of sample (c) was cho- 
sen to be comparable to the equivalent width threshold 
of 0.24 A from th e HST QSO absorption line key project 
(|Wevmann et al.|[l998l ). The last rows of the subtables 
in Table [7] show the mean dJ\f /dz for the 7 sightlines. 
We considered only absorbers with a\,lb,a^jN < 0.4. 
Note that the average values would have increased only 
by ~5% if we did not apply this cutoff. 

The average value of dN ' /dz in sample (c) (W > 0.24 
A) of 28 ± 4 for b < 40 kms -1 i s sligh tly smaller than 
the estim ates of IWevmann et aLI (|1998h at log(l + z) < 
0.15 and llmpev et alj (|1999f) at < z < 0.22, be- 
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TABLE 3 

Median, mean, dispersion of b at low z in a given column density interval for 

THE 7 QSO SAMPLE 



Column 
Density interval 



12.9,13.2 
13.2, 13.5' 
13.5, 13.8 
13.8, 14.1' 
14.1, 14.4' 
14.4, 14.7' 



median 


mean ±<r 


median 


mean ±<r 


median 


mean ±cr 


b > kms -1 


b < 40 kms -1 


b > 40 kms -1 


28.7 


31.5 ±13.9 


27.2 


26.7 ±6.5 


50.7 


54.7 ± 16.6 


34.0 


39.9 ±21.7 


27.0 


27.2 ±6.7 


55.4 


61.5 ± 21.1 


32.0 


36.3 ± 12.9 


28.7 


29.1 ±5.4 


54.1 


53.2 ±8.4 


34.0 


36.9 ± 14.1 


31.0 


30.7 ±4.9 


62.0 


57.8 ± 15.7 


37.0 


34.1 ±9.1 


27.7 


29.8 ±6.8 


45.7 


45.0 ± 1.2 


39.1 


35.2 ±8.9 


30.3 


30.5 ±8.0 


43.0 


43.2 ± 1.3 



Note. — Only data with a^/b, ctn/N < 0.4 are included in the various samples. 



TABLE 4 
Lyc* Line Density 



Siehtlinc 



Mn <W m /dz 
b < 40 kms -1 



A/hi dAfm/dz 
b < 150 kms -1 



(a) 13.2 < logAr m < 14.0 



H 1821+643 


7 


29 ± 11 


12 


50 ±15 


HE 0226-4110 


16 


40 ± 10 


22 


55 ± 12 


HS 0624+6907 


22 


67 ± 14 


31 


94 ±17 


PG 0953+415 


13 


64 ± 18 


17 


84 ±20 


PG 1116+215 


7 


56 ± 21 


12 


95 ±28 


PG 1259+593 


21 


68 ± 14 


35 


99 ± 17 


PKS 0405-123 


21 


51 ± 11 


35 


85 ± 14 


Mean 




54 ±5(15) 




80 ± 6(20) 


(b) 13.2 < logAT H i < 16.5 


H 1821+643 


9 


38 ± 13 


15 


63 ±16 


HE 0226-4110 


21 


52 ± 11 


31 


77 ±14 


HS 0624+6907 


25 


76 ± 15 


35 


106 ± 18 


PG 0953+415 


14 


69 ± 19 


18 


89 ±21 


PG 1116+215 


8 


64 ±23 


13 


103 ± 29 


PG 1259+593 


30 


85 ± 15 


43 


121 ± 19 


PKS 0405-123 


33 


80 ± 14 


47 


114 ± 17 


Mean 




66 ±6(17) 




96 ± 7(21) 


(c) 13.64 < logjVHi < 16.5 


H 1821+643 


6 


25 ± 10 


7 


29 ±11 


HE 0226-4110 


10 


25 ± 8 


16 


40 ±10 


HS 0624+6907 


10 


30 ± 10 


13 


40 ±11 


PG 0953+415 


6 


30 ± 12 


6 


30 ±11 


PG 1116+215 


2 


16 ± 11 


2 


16 ± 11 


PG 1259+593 


12 


34 ± 10 


18 


51 ± 13 


PKS 0405-123 


16 


39 ± 10 


21 


51 ± 12 


Mean 




28 ±4(8) 




37 ±4(13) 



Note. — Only data with ab/b, aj^/N < 0.4 are included 
in the various samples. Errors are from Poisson statistics. 
The number between parentheses in the row showing the 
mean corresponds to the standard deviation around the 
mean for the different lines of sight. 

cause the broader absorbers and uncertain absorbers are 
not included in our sa mple . Indeed, the estim ates of 
IWevmann et al.l (|1998l ) and llmpev et alj (119991) appear 
interm ediate between our two b samples. iPenton et alJ 
(|2004h found dN/dz = 25 ± (4, 5) for a sample with 
z < 0.069 which overlaps with our estimate within the la 
dispersion, which implies little or no evolution of dN / dz 
at z < 0.4 for the strong H I absorbers. For b < 150 
kms -1 , the value of dN/dz is about 1.3 times larger 
than for the NLAs. We find that the broad Lya lines 
with 40 < b < 150 kms -1 have dN(BLA)/dz « 9, im- 
plying that for the stronger lines of the Lya forest, the 
number of BLAs per unit redshift may be important. 
Comparison of samples (a) and (b) shows that the 



weak systems are far more frequent than the strong sys- 
tems. For these samples, we note that dhf j dz is system- 
atically smaller toward H 1821+643 and HE 0226-4110 
than dN/dz toward the other sightlines for both b sub- 
samples. dN/dz toward PG 1116+215 and PKS0405- 
123 is intermediate for the NLAs. The 3 other lines of 
sight have similar dN/dz. These trends do not appear 
related to the S /N of the data since the S /N is the highest 
toward H 1821+643 and comparatively low for HE 0226- 
4110. The redshift paths do not seem to explain all the 
differences. For example, while the redshift paths are 
comparable between HE 0226-4110 and PG 1259+593, 
Az is significantly smaller toward PG 1116+215. dN /dz 
is very similar for cither column density range toward 
the sightlines that cover small and large redshift paths, 
implying no redshift evolution of dN /dz between z > 0. 
and z < 0.4. Therefore, some of the observed variation 
in dN / dz must be cosmic variance between sightlines. 

We explore the effect of the S/N on the dN(BLA)/dz 
estimate in Table [5] by considering data with 
&b/b,aN /N < 0.4, 0.3, and 0.2. As expected, the spec- 
tra with the highest S/N are less affected by these cutoffs 
than the spectra with the lowest S/N. Decreasing the er- 
ror thresholds has an effect mostly on the weak systems 
(log iV ffl < 13.40). Toward PG 0953+415, the four BLAs 
have 13.2 < logiV H i < 13.40 and 0.23 < a b /b < 0.30, 
which explains why there is no BLA at the threshold 
a b /b < 0.2. Although the HS 0624+6207 spectrum has 
several weak BLAs, only if a b /b, <tn/N < 0.15 is applied 
do we observe a significant drop in dN (BLA) /dz by a fac- 
tor 2.5. On average, there is a decrease in dN(BLA)/dz 
by a factor ~1.4 between the cutoff at 0.4 and 0.2. If 
only BLAs with 40 < b < 100 kms -1 are considered, 
dN(BLA)/dz = 28,25,20 for a b /b,a N /N < 0.4,0.3,0.2, 
respecti vely. These e stimate s are in agreement with 
those of iRichter et~aT1 (|2006bD (see also gO}. The BLA 
density is about twice smaller than dN (NLA) /dz for 
samples (a) and (b). 

Finally, we summarize in Fig. [7] the frequency of the 
Lya absorbers with various 6-values in the low redshift 
universe. In this figure, dN/dz represents the mean of 
the number density of H I absorbers obtained toward 
each line of sight in four intervals of b ([0,20], [20,40], 
[40,60], [60,80] kms -1 ) using data with a b /b,a N /N < 
0.4 and logA+n > 13.2. The vertical error bars assume 
Poissonian errors. The mean dN /dz is shown at the 
mean &-value of each 6-interval sample. The horizontal 
bars show the fe-intcrvals delimited by the minimum and 
maximum ^-values of each ^-interval sample. The NLAs 
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TABLE 5 

Effect of the Signal-to-Noise on the Broad Lyo Density 
Estimate 



Sightline 


S/N 






cWni/rfz 








a b /b < 0.4 


a b /b < 0.3 


a b /b < 0.2 








a N /N < 0.4 


a N /N < 0.3 


a N /N < 0.2 


H 1821+643 


15 


-20 


25 ± 10 


25 ± 10 


25 ± 10 


HE 0226-4110 


5- 


11 


25 + 8 


20 + 7 


18 + 7 


HS 0624+6907 


8- 


12 


30 + 10 


30+ 10 


30+ 10 


PG 0953+415 


7- 


11 


20 + 10 


15 + 9 





PG 1116+215 


10 


15 


40 ± 18 


40+ 18 


32 ± 16 


PG 1259+593 


9 


17 


37 ± 10 


34 ± 10 


25 + 9 


PKS 0405-123 


5- 


10 


34 + 9 


24 + 8 


14 + 6 


Mean 






30 ± 4(7) 


27 + 4(8) 


21 ± 3(11) 



Note. — S/N is measured per resolution element. Sample with 40 < 
b < 150 kms~ 1 and 13.2 < log Nui < 16.5. Errors are from Poisson 
statistics. The number between parentheses in the row showing the mean 
corresponds to the standard deviation around the mean for the different 
lines of sight. 
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FlG. 7. — Mean number density per unit redshift with Poissonian 
errors (vertical bars) in four intervals of b ([0,20], [20,40], [40,60], 
[60,80] kms -1 ) using data with a b /b,a^ /N < 0.4 and logA^ni > 
13.2. The mean number density is shown at the mean fe-value 
of each 6-interval sample. The horizontal bars show the b-intervals 
delimited by the minimum and maximum b- values of each fe-interval 
sample. 

are mostly found between 20 and 40 kms -1 . Very nar- 
row Lya absorbers with b < 20 kms -1 are rare, and we 
note that in the [0,20] kms -1 interval, most of the ab- 
sorbers have 6 > 15 kms -1 . The paucity of very narrow 
absorbers with logiVm > 13.2 is real since the resolution 
of STIS E140M (6 inst - 4 kms" 1 ) allows to fully resolve 
absorbers with b > 5 kms -1 . The broad absorbers are 
more frequent in the b = [40,60] kms" range: there 
are about 3.5 times more absorbers in [40,60] kms" 1 
interval than in the b — [60, 80] kms" interval. Lya ab- 
sorbers are therefore more frequent for low column den- 
sities (log N(H I) < 14) and for b- values between 20 and 
40 kms 
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4. EVOLUTION OF THE DOPPLER PARAMETER 

4.1. Higher Redshift Samples 

Recently, the analysis of the Lya absorbers in the 
redshi ft range 0.5 < z < 2.0 from Jank necht et al.l 
(2006) has become available at the Centre de Donnees 
de Strasbourg (CDS). Their data that sample the mid-z 
IGM at 0.5 < z < 1.5 were obtained with HST/STIS 
E230M: PG 1634+706 (z = 0.534-1.295), PKS 0232-04 



(z = 0.876-1.419), PG 1630+377 (z = 0.875-1.451), 
PG 0117+213 {z = 0.875-1.475), HE 0515-4414 (z = 
0.874-1.475) and HS 0747+4259 [z = 0.760-1.443). The 
redshifts between parentheses indicate the redshift in- 
terval probed by the observations. The spectral reso- 
lution of E230M data (R ~ 30000) is lower than that 
of the low-z sample obtained with the E140M grating 
(R ~ 44000). The S/N of the mid-z sample is com- 
parable to the lowest S/N of the low-z sample, except 
for PG 1634+706, which has a S/N per resolution ele- 
ment of 5-40. The high redshift sample (z > 1.5) con- 
sists of data obtained with VLT/UVES (R ~ 40000): 
HE 0515-4414 (1.515-1.682), HE 0141-3932 (z = 1.518- 
1.784), HE 2225-2258 (z = 1.515-1.861), and HE0429- 
4901 (z = 1.662-1.910). The S/N is typically higher than 
30-40, except for HE 0429-4901, which is about 15. The 
spectrum of HS 0747+4259 (z = 1.562-1.866) was also 
obtained with Keck/HIRES (R ~ 50000) with a S/N in 
the range of 6-24. 

For the high redshift sample, we also consider the 
spectra of QSOs at z > 1.5 presented by Kim et al. 
(2002a). Their profile fitting results are available at the 
CDS. These data were obtained with the VLT/UVES 
at resolution of R ~ 45000 and typical S/N~ 40-50. 
The QSOs considered are: HE 0515-4414 (z = 1.53- 
1.69), Q 1101-264 (z = 1.66-2.08), J 2233-606 (z = 1.80- 
2.20), HE 1122-1648 (z = 1.88-2.37), HE 2217-2818 (z = 

I. 89-2.37), HE 1347-2457 (z = 2.09-2.57), Q 0302-003 
(z = 2.96-3.24), Q 0055-269 {z = 2 99-3.60). Note 
that t he HE 0515-4414 VLT spectrum in lJanknecht et all 
( 2006) has a different wavelength coverage and exposure 
ti me than the HE 05 15-4414 VLT spectrum presented 
bv iKim et al.l (|2002ah . We finally consider the H I pa- 
rameters measured toward HS 194 6+7658 [z = 2.43- 
3.05) by iKirkman k Tvtlerl (|1997l ). The spectrum of 
HS 1946+7658 was obtained with Keck/HIRES. The S/N 
per resolution element varies from about 30 to 200 and 
the spectral resolution is similar to UVES. The spectral 
resolution of the high redshift sample is similar to the low 
redshift sample, but the S/N is generally much higher in 
the z > 1.5 sample than in the low or mid z samples 

4.2. Conditions for Comparing Various Samples 

The definition of a BLA that was followed in recent 
papers presented by our group is an absorber that can 
be fitted with a single Gaussian component with b > 40 
kms" 1 and for which the reduced-% 2 does not improve 
statistically by adding more components in the model. 
Low S /N data can have, however, treacherous effects that 
can mask a BLA or confuse narrow multi-absorbers with 
a BLA (see Figure 2 in Richter et al. 2006b). To over- 
co me the effects of nois e, the sample of BLAs defined 
bv iRichter et al.l ()2006bf ) has also to satisfy the follow- 
ing rules: (i) the line does not show any asymmetry, (ii) 
the line is not blended, (iii) there is no evidence of mul- 
tiple components in the profile; (iv) the S/N must be 
high enough (i.e. log(iV H i/6) > log [3 x 10 12 /(S/N)] > 

II. 3, where Nm and b are in cm" 2 and kms -1 , re- 
spectively). For the low redshift sample presented here, 
the BLAs strictly follow criteria (iv) when the condi- 
tion aft/b, on /N < 0.4 is set.. Following these rules, 
IRichter et all (|2006bD found dN '(BLA) /dz = 22 ± 5 for 
their secure detections (53 for the entire candidate sam- 
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pie). Within la, our BL A number de n sity es timate 
overlaps with the result of iRichter et alJ (|2006bl ) when 
the cutoff Ob/b,aN /N < 0.4 is applied to our sample. 
Therefore, this shows that by applying an error cutoff 
without scrutinizing each profile for the conditions listed 
above, we find a similar average cW(BLA) jdz. This 
is crucial because for a comparison with other samples 
at higher z, we can not examine eac h absorber indi- 
vidual ly. For the sample presented by Uanknecht et all 
(2006) , no spectra or fi t s are shown. iKim et al.l (120023) 
and lKirkman fc Tvtlerl (|1997l ) present their spectra and 
fits, but over too broad wavelength ranges to study in 
detail the conditions listed above. Furthermore, as z 
increases absorbers are more often blended due to the 
higher redshift line density, so we can not reject the 
blended systems, otherwise we would introduce a strong 
bias in the comparison. At high, mid, and low redshift, 
the x 2 °f the profile fit governs the number of compo- 
nents allowed in the model of an absorber; therefore a 
similar methodology was applied in each sample, allow- 
ing a direct comparison of the various measurements. We 
apply the same cutoff crjy/N < 0.4 to the z > 0.5 
redshift samples to remove the uncertain profile fit re- 
sults in a similar manner in each sample. Such cutoff, 
however, introduces a systematic effect: more BLAs will 
be rejected in low S/N spectra (low and some mid z data) 
than in high S/N spectra (high z data). But such sys- 
tematics should underestimate the number of BLAs in 
low S/N data, and therefore this should strengthen the 
differences observed at low z compared to the higher z 
samples. 

For our comparison, we also consider only absorbers 
with log A/hi > 13.2, the completeness level of the low 
redshift sample, which also corresponds to about the 
completeness of the lowest S/N spectra of the mid-z sam- 
ple. If the S/N of the high redshift spectra is solely 
considered, absorbers with log A/hi > 13.2 would be far 
above the completeness level of the high redshift sam- 
ple, which is logA Hi ^ 12.5 for the data presen ted by 
IKim et all (j2002al ) and IKirkman fc Tvtlerl (|1997l ). How- 
ever, line blending and blanketing re duce the complete- 
ness t hreshold, espec ially at z > 2.5 (Kirkma n fc Tvtlerl 
[19971 : ILu et al.lfl99ll . At z ~ 4, ILu et al.l (119961 ) show 
using simulated spectra that line blanketing was not 
imp ortant as long as l og Afa i > 13.5, while at z ~ 
2.7, IKirkman fc Tvtlerl (|1997l ) show following a similar 
methodology that their sample is likely to be complete 
at logA^Hi > 12.8-13.0. Since broad and shallow ab- 
sorbers are more uncertain, considering only absorbers 
with log A/hi > 13-2 also reduces the problem of creating 
a higher proportion of wide lines than is really present in 
the intrinsic distribution. 

For our comparison, we only consider BLAs with 40 < 
b < 100 kms -1 . The choice of b < 100 kms -1 reduces 
the incompleteness at the high b end of the low red- 
shift sample. Furthermore, at z ~ 4 (which is at higher 
redshift than any absorbers considered here), ILu et all 
(1996) argue that absorbers with b > 100 kms are 
caused essentially by heavily blended forest lines since 
they are systematically found in the high line density re- 
gion of the spectrum. This eff ect should diminish g reatly 
as z decreases. At z ~ 2.7, IKirkman fc Tvtlerl l| 19971 ) 
have produced simulated spectra to better understand 



the intrinsic properties of their observational data. They 
observe a tail at high b in both the simulated and ob- 
served b distributions. In the simulated spectra, the tail 
at b > 80 kms -1 is only due to line blending because 
their simulation did not allow such broad lines. They 
find more lines at b > 80 kms -1 in the observed distri- 
bution (3.6% compared to 0.5% in the simulated spec- 
tra), suggesting that many of these absorbers could be 
intrinsically broad. Refined simulations, so that simu- 
lated observations match exactly the real observations, 
would be needed to be e ntirely conclusive on this point 
(|Kirkman fc Tvtlerl [i~997t) . Yet, at high z, if BLAs with 
b > 80 kms -1 that are in fact blended narrow lines 
were frequent, this effect should be more important as 
z increases. We will see that this effect, if present, is 
not statistically significant (see next two sections). By 
considering only absorbers with 40 < b < 100 kms , 
log A/hi > 13.2, and (Tb/b, ctn/N < 0.4, we reduce signifi- 
cantly the risk of including spurious broad absorbers. We 
also note that, in any samples considered, the majority 
of BLAs are actually found in the 6-range 40 < b < 60 
kms - , not in the 6-range b > 80 kms - . 

In Fig. [5J we show the 6- A/hi distribution for three 
redshift intervals, z < 0.5 (top panel), 0.5 < z < 1.5 
(middle panel), and 1.5 < z < 3.6 (bottom panel). Only 
absorbers that satisfy the conditions < b < 100 kms -1 , 
log A/hi > 13.2, (Th/6, ctjv /N < 0.4 are shown in the fig- 
ure. At z > 0.5, there are many weak systems with 
b < 40 kms -1 , which darkens part of the middle and 
bottom panels. At z < 0.5, there are very few absorbers 
with b > 40 kms -1 and log A^hi > 14, and none of them 
have b > 50 kms -1 . This contrasts remarkably with 
what is observed at higher redshifts: many absorbers 
at z > 0.5 have b > 40 kms -1 and log Am > 14.0. 
This effect must be due to strong, saturated Lya lines 
for which the errors for b and logA^i are far too op- 
timistic. At z < 0.5, saturation can be dealt with be- 
cause higher Lyman series lines are systematically used 
in the profile fitting, reducing the possibility of non- 
uniqueness in the profile fitting solution and of finding 
unrealistic large b- values for strong lines. At z > 0.5, 
Lya is g enerally th e only t ransition available, although 
we note IKim et ail (|2002aD analyzed the Ly/3 forest at 
z < 2.5 with a small sample and their results suggested 
that line blending and saturation are not an important 
issue. Yet, Fig. [5] shows a system with b = 97.0 ± 3.8 
kms -1 and logA^i = 16.44 ± 0.14. This absorber is ob- 
served in the spectrum of J 2233-606 at z — 1.869, and 
the parameter s that a re presented here were estimated 
by IKim et all (l2002af). The sp ectrum of J 2233-606 is 
shown in iCristiani fc D'Odoricol (pOOl and at the wave- 
length correspondi ng to this absorber , there is a strong, 
saturated line th at IKim et al.l (l2002al) model w ith a sin- 
gle component. iD'Odorico fc Petitjeanl (j2001h actually 
found that this absorber has ~11— 16 components using 
higher Lyman series lines and metal lines. While this ex- 
ample is extreme, most of the absorbers with log A/hi > 
14 and b > 40 kms -1 are lik ely to be strong sa t urate d 
Lya lines at z > 0.5. Fig. 1 in IKirkman fc Tvtlerl (|1997l) . 
which shows the whole HS 1946+7658 spectrum with the 
values of log A/hi, 6, corroborates this conclusion. The 
strong, saturated Lya lines are unlikely to probe a single 
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Fig. 8. — Left panel: Line width versus H I column density for the Lya absorbers for 13.2 < logAfjii < 16.5 and < b < 100 kms — 1 
in different redshift intervals indicated in the upper-right corner. Only data with less than 40% errors in b and TVpji are considered. The 
symbols have the following meaning: data represe nted by circles we r e esti mated from the low redshift sample (see Table[l]for the references); 
data represented by stars are from the sample of Janknccht et a . (2006); data represe nted by triangles are from the sample of Kim et al. 
(2002a); data represented by squares are from the sample of Kirkman & Tvtlcr (1997). Right panel: Same as the left-hand side, but only 
for absorbers with 13.2 < log JV H i < 14.0 



broad line, but are likely composed of several unresolved 
components. Therefore, for our comparison, we will con- 
sider absorbers with 13.2 < logiVni < 14.0. This sample 
is summarized in the right-hand side of Fig. [H] 

As we have just illustrated, the Voigt profile fitting 
method adopted for deriving z, 6, A/hi may not be unique, 
especially for highly blended regions, strong lines, and 
low S/N spectra. We believe that some of this effect is 
reduced by considering absorbers with 13.2 < logA^ni < 
14.0, b < 100 kms -1 , and a b /b,a N /N < 0.4. At low 
z, line blending and strong lines are not a major issue 
because often higher Lyman series lines can be used. 
Futhcrmore, several persons analyzed these data inde- 
pendently, and other methods (curve-of-growth and/or 
optical depth method) were also used for the low-z sam- 
ple, yielding consistent results (see the references given in 
Table [1]), although the Appendix highlights some differ- 
ences between various groups. Even in this latter case, 
there is roughly a la consistency between the various 
results f or absorbers with 1 3.2 < . logAni < 14.0. At 
z ~ 2.7, iKirkman fc Tvtleri (|1997l) show that profile fit- 
ting may not provide a unique solution even in very high 
S/N data. Yet, their comparison with simulated spectra 
show that the intrinsic distribution of b is very similar 
to the observed b distribution. Therefore, in a statistical 



sense, the observed b distribution at high z can be used 
for our comparison. We also note that three different re- 
search groups have worked on the absorbers at z > 1.5 
presented here, with the spectra obtained from different 
telescopes with different S/N, and we did not find any 
major differences between their results, at least within 
the conditions listed above. 

Continuum placement may also be a problem for 
finding or defining a BLA. At low z, BLAs may be 
confused with continuum undulations. With our con- 
strai nts, most of these i ll-defined BLAs are rejected (see 
also iRichter et al.1 l2006bf ) . We believe that this con- 
clusion should apply to the mid-z sample, but since 
Janknccht et al.l (|2006[ ) did not present any of the spec- 
tra they analyze, we treat this sample with more caution 
(see also below). At z > 1.5, the continuum is more diffi- 
cult because each order must be fitted with a high order 
Legendre polynomial, and during the Voigt profile fitting 
process, the continuum level is often adjusted to give an 
optimum fit between the data and the calculated Voigt 
model. Furthermore, as z increases the continuum place- 
ment becomes more difficult because the line density in- 
creases, and therefore the line blending and blanketing 
effects increase. High order polynomials and adjusted 
continua can therefore mask broad and shallow BLAs. 
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However, such effects may be counterbalanced by a pos- 
sible spurious increase of BLAs caused by the same line 
blending and blanketing effects (see above). The high 
S/N of the data at high redshifts and considering only 
the absorbers with 13.2 < logiV H i < 14.0 and b < 100 
kms -1 reduce greatly the risk of losing many BLAs be- 
cause of the line blending and blanketing at high z. The 
problem of not finding BLAs is potentially more impor- 
tant in the highest z spectra, but we do not notice signif- 
icant differences between z ~ 2 and z > 3. We also note 
above that simulated spectra show that line blanketing 
was not important at logTVni > 13.2. Therefore, for the 
samples considered here using the limits log Nm > 13.2 
and b < 100 kms - , they must almost be complete. Ul- 
timately, one would like to produce simulated spectra of 
various z with realistic inputs to fully consider the im- 
pacts of S/N, continuum placement, and line blending 
and blanketing on finding BLAs and deriving their in- 
trinsic properties. 

The redshift range z = 0.5-1.5 is more problematic 
than the other redshift ranges because most of the data 
have low S/N, all have lower resolution than the higher 
or lower r edshift spectra, and no s pectra or fits were pre- 
sented by iJanknecht et al.l (|2006f ) making it more diffi- 
cult to assess some of the issues discussed above. Ly/3 
lines are only available for about half the wavelength 
coverage of this redshift range and lower spectral resolu- 
tion and S/N of these data further worsen the problems 
of saturation and non-uniqueness in the profile fitting 
results. Setting the conditions < b < 100 kms , 
13.2 < logJVki < 14.0, and a b /b,a N /N < 0.4 make 
this sample stronger, but we nonetheless treat the re- 
sults from this sample with caution. We will consider 
two sub-sam ples for the estim a te of the redshift density 
based on the IJanknecht et al.l ([2006) sample (see §4.4|) : 
one with all their data, and one with only their highest 
S/N data (i.e. PG 1634+706, HE 0515-4144, HE 0141- 
3932, HE 2225-2258). Among the high S/N data in the 
mid-z sample are E230M spectra of PG 1634+706 and 
HE 0515-4144. 

For the reasons aforementioned, with the conditions 
< b < 100 kms" 1 , 13.2 < logA H i < 14.0, and 
<7b/b,<7N /N < 0.4, we greatly reduce some pitfalls in 
comparing data analyzed by various groups, with dif- 
ferent S/N, obtained with various instruments, with dif- 
ferent line blending issues. A BLA is therefore defined 
here as an absorber that is fitted with a single Gaussian 
for which the \ 2 does not drop significantly by adding 
more components (this condition was adopted by the var- 
ious groups who analyzed the data used here) and that 
has 40 < b < 100 kms" 1 , 13.2 < logiV H i < 14.0, and 
0&/6, <7n/N < 0.4. We believe that we are statistically 
comparing the intrinsic properties of the broadening of 
H I absorbers, although we note that simulated spec- 
tra probing various z, with realistic inputs may be the 
only way to fully understand some of the effects discussed 
above and how these effects balance each other at various 
z. 

4.3. Comparison of the Distributions of b at Low and 

Higher z 

Table [5] summarizes the b- value median, mean, disper- 
sion, and fraction of BLAs (i.e. systems with 40 < b < 



100 kms 1 ). In this table, we consider two sub-samples 
with < b < 100 kms - : (a) the entire sample, (b) 
the higher quality sample with 13.2 < logiVni < 14.0 
and errors in b and ./Vhi less than 40%. We note that 
the number of systems in e ach sam ple is roughly similar, 
except for the iKim et ail (|2002a| ) sample, which is no- 
ticeably larger, and the 0.5 < z < 1.0 sample, which is 
somewhat smaller. In sample (a) there appears to be an 
increase in the median, mean, and the fraction of BLAs 
as z decreases. The same trend is observed in sample (b), 
but the differences are better revealed, where the median 
and mean are always larger by 15-30% in the low red- 
shift sample (z < 0.5) than in the high redshift sample 
(z > 1.5). For the sample (b), the fraction of BLAs at 
z < 0.4 is 1.9-2.4 times larger than the fraction of BLAs 
at z > 1.5. Therefore, b is larger on average at z < 1.0 
than at z > 1.5, and the fraction of BLAs is larger at 
z < 1.5 than at z > 1.5. 

In Figs.l9TfT2"lwe compare the distributions of 6(H I). In 
all comparisons we restrict the observations to the higher 
quality measurements with less than 40% errors in b and 
N m . 

In Fig. [9] we compare the normalized number dis- 
tribution of 6(H I) for the low redshift sample with 
the normalized numbe r distribution of 6(H I) for the 
IJanknecht et al.l (|2006h sample, which is subdivided in 3 
z sub-samples. The distributions peak at about b ~ 20- 
30 kms -1 . Each distribution shows a tail at higher 
b. The tail in the 1.5 < z < 2.0 distribution is al- 
ways weaker than in the lower z-interval distributions, 
especially for 40 < b < 55 kms -1 . The tail in the 
0.5 < z < 1.0 distribution is generally stronger than 
any other ^-interval distribution presented in this figure, 
which is possibly due to a combination of blending effects, 
low S/N spectra, and lower resolution spectra in the red- 
shift range 0.5 < z < 1.0. We note that data probing 
the redshift ranges 0.5 < z < 1.0 and 1.0 < z < 1.5 were 
obtained with the STIS E230M grating, but the sensi- 
tivity of E230M is lower for Lya in the redshift range 
0-5 < z < 0.9 than 0.9 < z < 1.4. A Kolmogorov- 
Smirnov test does not reveal a significant difference be- 
tween the samples z < 0.5 and 0.5 < z < 1.5 (the maxi- 
mum deviation between the two cumulative distributions 
is D = 0.118 with a level of significance P — 0.104), but 
suggests a difference between the samples z < 0.5 and 
1.5 < z < 2.0. Therefore, BLAs appear more frequent at 
z < 1.5 than at z > 1.5. 

In Fig. 1 101 we compare the normalized number dis- 
tribution of 6(H I) for the low redshift sample with the 
normalized number distribution of 6(H I) for the high 
redshift sample (1.5 < z < 3.6) from Kim et al. (2002a). 
The left panel shows the entire sample at 1.5 < z < 3.6, 
while the right panel shows the distribution of b in var- 
ious redshift intervals for the high z sample. The peak 
of the high redshift normalized distribution in Fig. [10] 
is not only shifted by ~ —5 kms -1 with respect to the 
peak of the low redshift distribution, but also the width 
of the distribution is smaller for the high redshift sam- 
ple. Moreover, the low redshift sample shows a tail of 
high 6(H I) absorbers that is much weaker in the high 
redshift sample, showing that BLAs are more frequent at 
z < 0.5 than at z > 1.5. The right-hand side of Fig. [10] 
verifies this conclusion for the various redshift intervals, 
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TABLE 6 
Evolution of b 





z < 0.5 
(this paper) 


0.5 < 2 < 1.0 1.0<2<1.5 
(J06) (J06) 


1.5 < 2 < 2.0 
(J06) 


1.5 < 2 < 3.6 
(K02) 


2.4 < 2 < 3.1 
(KT97) 








< b < 100 kms" 


l 






b median (m to t) 
b mean ±<r 
m(b > 40)/m to t 


30.5 (336) 
34.0 ± 16.9 
0.277 


27.9 (263) 28.3 (588) 
32.9 ±17.8 31.8 ±17.0 
0.274 0.240 


28.0 (450) 
31.4 ± 14.5 
0.204 


24.8 (2305) 
27.2 ± 13.5 
0.131 


27.1 (452) 

31.2 ± 16.2 
0.210 


< b < 100 kms -1 and < (a b /b,a N /N) < 0.4 and 13.2 < logiYm < 14.0 


b median (mtot) 
b mean ±cr 
m(b > 40)/m tot 


32.7 (162) 
37.4 ± 16.1 
0.321 


31.5 (82) 30.8 (209) 
36.0 ± 18.3 35.5 ±17.3 
0.317 0.272 


28.2 (154) 
31.7 ± 12.9 
0.169 


25.5 (509) 
29.2 ± 11.9 
0.132 


27.1 (123) 
32.5 ± 16.0 
0.179 



Note. — Median and mean values of b are listed for the different redshift intervals for the entire samples with 
< b < 100 kms -1 and the higher quality samples restricted to absorbers with 13.2 < logiVHi < 14.0. mtot is the 
total number of absorbers in the sample. m(b > 40) is the total number of absorbers with 40 < b < 100 kms" 1 in the 
sam ple. The values o f er lis ted are the standard dev iation ar ound the mean. Source o f profile fitting measurements: 
J06: IJanknecht et al.l J2006I) : K02: IKim et al.l <l2002aD ; KT97: IKirkman fc Tvtlerl U997D . 
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Fig. 9. — Comparison of the normalized distributions of b(H I) 
for the low 2 sample (bl ack histogram) and for t he higher 2 sample 
(colored histograms) of Janknccht et al. (2006) for the weak ab- 
sorbers 13.2 < logA^i < 14.0. The sample in the redshift range 
2 ~ 0.5—1.5 has a spectral resolution of about 10 kms -1 , while the 
other samples have about 7-8 kms -1 . Only data with less than 
40% errors in b and JVhi are considered. 

i.e. BLAs appear more frequent at z < 0.5 than for any 
redshift intervals at z > 1.5. There is, however, no ma- 
jor difference between the various redshift intervals at 
z > 1.5. We also note that the b distribution at z > 1.5 
shows little effect of line blending as z increases since 
there is scant evidence of a larger of fraction of BLAs at 
z > 3 than at 1.5 < z < 2. 

We compare in Fig.[Tl]our low z sample to the high red- 
shift absorbers observed toward the QSO HS 1946+7658 
(jKirkman fc Tvtlerl [l997h . The median and mean are 
somewhat interm ediate between our sample and the 
IKim et al.l (|2002al ) sample (see Table [6]). We slightly ad- 
justed the redshift of the low z sample so that the low 
and high redshift samples have exactly the same num- 
ber of systems (m — 130); hence the two distributions 
can be directly compared. The peak of the high red- 
shift distribution in Fig. [10] is again shifted by ~ — 5 
kms -1 with respect to the peak of the low z distribu- 



tion. A higher number of systems with 40 < b < 70 
kms -1 is found in the low redshift sample. At b > 70 
kms -1 the high redshift sample appears to have a larger 
number of systems . How ever, we note that the data of 
IKirkman fc Tvtlerl (|1997f ) have S/N up to 200, and as 
we discussed in our sample is not complete for b > 80 
kms -1 . A Kolmogorov-Smirnov test yields a maximum 
deviation between the two cumulative distributions of 
the low and high redshift samples D — 0.210 with a level 
of significance P = 0.006; the null hypothesis of no dif- 
ference between the two datasets is therefore rejected. 

So far, we have ignored the effects of evolution of col- 
umn density in the absorbers in our c omparison. Accord - 
ing to the numerical simulation of iDave et alJ (l999), 
the dynamical state of an absorber depends mainly on 
its overdensity, so that a column density range traces 
absorbers of progressively higher overdensity as the uni- 
verse expands. Therefore, according to Eq. [21 a low red- 
shift H I absorber is physically analogous to a high red- 
shift H I absorber not with the same column density but 
to an absorbe r with column de nsity l0 - 4x( ^ hi e h-2l ° w )/ - 7 
times higher (jDave eFaDEMi). With the average red- 
shifts of 2 . 6 and 0.2 of the high redshift sample of 
IKim et ail (|2002af) and our sample, this corresponds to 
a column density roughly 25 times higher, i.e. a sample 
with 12.5 < logAni < 14.1 at low z is physically analo- 
gous to the high redshift sample with 13.9 < logAni < 
15.5. At such high column density, BLAs are likely to 
trace in large part narrower, strong absorbers that are 
blended together (see §4.2|) . We, nonetheless, make this 
comparison in Fig. 1121 where both samples have the same 
number of systems. The high redshift sample has again a 
larger fraction of systems in the range b — [20, 40] kms - 
than the low redshift sample. The fraction of BLAs at 
low redshift is again larger than the fraction of BLAs 
at high redshift with b > 50 kms -1 . This is remark- 
able since the low redshift sample is far from complete 
for logA^i < 13.2 and the high redshift sample likely 
overestimates the number of true BLAs (i.e. BLAs that 
are not blends of NLAs) . We again apply a Kolmogorov- 
Smirnov test on the two cumulative distributions of the 
low and high redshift samples and find D — 0.130 with a 
level of significance P = 0.029; again the null hypothesis 
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Fig. 10. — Left panel: Comparison of the norm alized distr i butions of f>(H I) for the low z sample of data points (black histogram) and 
for the high redshift sample (re d histogram) from Kim ct al. (2002a) for the weak absorbers 13.2 < log^Hi < 14.0. Right panel: Same as 
the left panel but the sample of Kim et al. ( 2002a]) is separated in various redshift intervals. Each sample has a similar spectral resolution 
of about 7-8 kms . Only data with less than 40% errors in b and JVhi are considered for the different samples. 
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Fig. 11. — Comparison of the distributions of 6(H I) for the low 
z STIS E140M sample (bl ack histogram) and for th e high z Keck 
sample (red histogram) of Kirkman & Tytler ( 1997) for the weak 
absorbers 13.2 < log TVhi < 14.0. The number of absorbers in each 
sample is the same: m = 130. Only data with less than 40% errors 
in b and TVjji are considered for both samples. 

of no difference between the two datascts is rejected. 

4.4. Evolution of dJ\f/dz as a Function of b 

In Fig. [H we show dAf /dz (0 < b < 100 kms" 1 ), 
dAf (NLA) /dz (b < 40 kms" 1 ), dAf(BLA)/dz (40 < b < 
100 kms" 1 ), and {dAf (BLA) / dz) /(dAf '(NLA) i/ dz) as a 
function of the redshift. Only systems with 13.2 < 
logA^Hi < 14.0 and ab/b,aN/N < 0.4 are considered 
in this figure. dAf/dz was estimated for each sightlinc 
over the redshift interval available along a given sight- 
line. The redshift at which dAf /dz is plotted in Fig. [T3] 
corresponds to the mean redshift interval in a given sight- 
line. The vertical bars are Poissonian errors, while the 
horizontal bars represent the standard deviation around 



Fig. 12. — Comparison of the distribution of 6(H I) for the low z 
STIS E140M sample of data points (black histogram) and for the 
high redshift sample (red histogram). The samples cover different 
ranges in logiVm. Higher column density systems at high redshift 
are expected to be physically analogous to the low density systems 
at low redshift in an expanding universe. But note that in the high 
redshift sample, BLAs are often likely to b e blends of narrowe r lines 
(see §4.2). The high redshift data are from lKim et al.l (I2002all . The 
number of absorbers in each sample is the same: m = 230. Only 
data with less than 40% errors in b and TVjji are considered in both 
samples. 

the mean of the ob served redshifts of the L ya absorbers. 
For the sample of Jankncc ht et al.l (|2006( ), several lines 
of sight have a redshift path larger than 0.6 and we di- 
vided those in two redshift intervals. The left diagram 
in Fig. [T3l includes all the sightlines available in the var- 
ious samples, while the right diagram includes all the 
sightli nes available in the low z sample and high z sam - 
ples of lKim et all (|2002af) and lKirkman fc Tvtlerl (11997ft. 
but on ly the highest quality data in the I Janknecht et all 
(2006) sample (see JO]). 
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Fig. 13.— The Lya redshift density of the entire sample (0 < b < 100 kms -1 ), the NLAs (6 < 40 kms -1 ), the BLAs (40 < b < 100 
kms" ), and the ratio of the Lya redshift density fraction of the BLAs to the NLAs (bottom panel). We only consider absorbers with 
13.2 < log JVhi < 14.0 a nd with errors less than 40% in b and Nm. In the top two panels, the solid line shows dj\f /dz <x (1 + z) 1,42 adopted 
from lKim et al. (2002a|). The solid gray crosses in the bottom panels are the mean of (eiA/*(BLA) /dz)/(dj\f (NLA) /dz) where the horizontal 
bars indicate the redshift range over which the mean is estimated. The vertical bars are Poissonian errors. The symbols have the following 
meaning: data represented by ci rcles are estimated from the low redshift sample (see Table [T] for the references); data represented by star s 
are estimated from the sample of Janknecht et al. (2006); data re presented by triangles are e stimated from the sample of [Kim et al. ( 2002a); 
data represented by the square are estimated from the sample of Kirkman & Tvtler (1997). The left panel considers all t he data with the 
condit ions listed above from these various samples, while in the right panel we remove the lowest quality measurements of Janknec ht et al] 
i2006h (see §4.2). 



The top two panels of Fig. [T3] show the usual num- 
ber density evolution with little or no evolution between 
rcdshifts and ~1.0-1.6, and an evolution of dM/dz 
with z at higher redshift (see for example Kim et al. 
2002a). At high redshift, dM /dz decreases with decreas- 
ing z (dM/dz oc (1 + z) 7 , 7 > 0) according to the ex- 
pansion of the universe, wh ich forces any ini tial baryon 
overdensity to thin out (e.g.. lDave et al.|[l999t and refer- 
ences therein). The expansion also results in a decrease 
of recombinations of the free electrons with the protons 
and thus in an additional decline in dM/dz. The break 
in dM /dz near z ~ 1.6 (log(l + z) ~ 0.4) is believed to be 
primarily caused by the drop in the UV b ackground be- 
cause of the declining quasar population (|Theuns et al.l 
H998U D ave et aLlll999f ). In the top panels, the solid line 
shows dM /dz = (dM /dz ) n(l + z ) 7 , where 7 = 1.42 was 
adopted from lKim et al . (2002a) who find this value for 
the weak absorbers with 13.1 < logiVm < 14.0. The 
value of (dM /dz)o was adjusted to match the data. For 
the entire sample (0 < b < 100 kms -1 ) or the sample 
restricted to b < 40 kms -1 , this line represents well the 
evolution of dM /dz at z > 1, but we note that when the 
best quality data are considered (right-hand side panel), 
the break in the evolution appears to occur at z ~ 1.6. 



The second panels from the bottom in Fig. [13] show 
that dM(BLA)/dz is generally similar or higher at low 
and mid redshift than between redshifts 1.6 and 2. At 
z > 2, dM (BLA) /dz increases and is larger at z > 2.5 
than dM(BLA)/dz at z < 1.5. The Hubble expansion 
must be the primary driver in the evolution of these 
structures at high redshift, and despite the expansion 
the BL A number density is comparable to d/V(BLA) / dz 
at high z, which contrasts remarkably from the evolu- 
tion of dM (NLA) /dz from low to high z. This differ- 
ence is even more striking when we consider the ratio 
K = (dM (BLA)/ dz) / (dM (NLA) /dz) (bottom panels), 
which is higher in the low redshift sample than in the 
high redshift sample. On the left-hand side diagram, at 
mid redshift there is a very large scatter in the redshift 
range 0.9 < z < 1.2. This scatter clears up when only the 
best quality data of the mid-z sample are considered (see 
right-hand side diagram) . The distribution of 1Z between 
rcdshifts ^1.6 to 3.5 is actually nearly flat, showing that 
the slope controlling the evolution of d/V(NLA) /dz and 
dM '(BLA) / dz at z > 1.6 must be about the same, further 
strengthening that NLAs and BLAs must follow a similar 
evolution dictated by the expansion of the universe. We 
note that TZ(z > 2.5) is slightly larger than TZ(z ~ 2), 
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possibly indicating some spurious increase of BLAs pos- 
sibly due to line blending and blanketing effects that are 
more serious at z > 2.5 than at z ~ 2. However, in 
the view of the large scatter in the various samples, this 
result does not appear statistically significant. The grey 
crosses in the bottom panels show the average 1Z over the 
redshift range depicted by the horizontal grey bar (the 
vertical grey bar assumes Poissonian errors). At z < 1.6, 
1Z is larger than at z > 1.6 as both the individual 1Z and 
average 1Z show. 7Z is a factor 3.2 highe r at z < 0.4 than 
at z > 1.6. When only the best data of iJanknecht et all 
(2006) are considered, TZ at mid-redshift (0.6 < z < 1.5) 
appears similar to the value observed at low redshift. 

If we had considered systems with b < 80 kins" 1 , the 
same conclusions would be drawn. If we had considered 
a cutoff for the NLAs of b — 50 or 60 kms -1 instead 
of 40 kms -1 , similar conclusions would also be drawn. 
If absorbers with 13.2 < log iVni < 16.5 are considered, 
we find that 7Z(z < 0.5) ~ 3TZ(z > 1.6), despite the fact 
many BLAs at high z may actually be blends of narrower 
lines (see JO]). 

4.5. Implications 

We find that (1) the distribution of b for the BLAs 
has a distinctly more prominent high velocity tail at low 
and mid z than at high z; (2) the median and mean b- 
values at low and mid z are systematically higher than at 
high z; and (3) the ratio (dN (BLA) /dz) / (dJ\f (NLA) / dz) 
at low and mid z is a factor ~3 higher than at high z. 
These conclusions hold for a division b = 40, 50, or 60 
kms -1 between NLAs and BLAs. For the reasons dis- 
cussed in §4.2[ we believe that these are intrinsic proper- 
ties of the evolution of the physical state of the gas that 
are not caused by comparing data with different S/N and 
line blending and blanketing effects. However, simulated 
spectra probing various z, with realistic inputs would cer- 
tainly help to unravel how exactly some of these issues 
(e.g., continuum placement, line blending effect, differ- 
ent S/N ratios) balance each other. Our results strongly 
suggest that if the broadening is mostly thermal, a larger 
fraction of the low z universe is hotter than the high z 
universe, and if the broadening is mostly non-thermal, 
the low z universe is more kinematically disturbed than 
the high z universe. It is likely that both possibilities are 
true. 

5. THE DIFFERENTIAL COLUMN DENSITY 
DISTRIBUTION FUNCTION 

The differential column density distribution f(Nm) is 
defined such that f(Nm,X)dXdNm is the number of ab- 
sorption systems with column density between Ajji and 
Ahi + d/Vffl and redshift path between X and X + dX 
(e.g.. lTvtledli987h . 

Tfi 

f(N m )dN m dX = AN j, Ax dNmdX , (5) 

where m is the observed number of absorption systems 
in a column density range AAhi centered on Am ob- 
tained from our sample of 7 QSOs with a total ab- 
sorption distance coverage SAX = 2.404 (see $2|). 
Empirically, it has been shown that at low and high 
redsh ift, /(Ahi) i s well fitted by a single power law 
(e.g., lTvtlerl[l987t iPetitiean et ailfl99l iHu et alJH995t 
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Fig. 14. — The differential density distribution (f(Nui)) without 
an incompleteness correction is plotted against log TVjji • Only H I 
absorbers with fe(H I) < 40 kms -1 and errors less than 40% in 
b and TVjji are considered. The data are binned only for purpose 
of presentation. The solid line is a maximum-likelihood fit to the 
data (f(N m ) = C K iN~f) with the slope = 1.76 for H I column 
density systems between 13.2 and 16.5 dex. The dotted and dashed 
lines are maximum-likelihood fit to the data for 13.2 < logAfjji < 
14.4 and 14.4 < logJV H l < 16.5, respectively (see Table [7] and §S\ 
for more details). 

Lu et all H996L lKimgt~afl l2002at iPenton et all I2000L 
200i iDave fc Trindl27)0iTk 

f(N m )dN m dX = C m N-fdN m dX . (6) 

In Table [7J we show the results from the maximum- 
likelihood estimate of the parameters for the slope /3 
and the normalization constant where Chi = Wtot(l ~ 
/3)/[A^(l - (A r min/A ma ,) 1 -' 3 )] where m tot is the to- 
tal number of absorbers in the column density range 
[A m i n ,A max ]. We separate our sample into absorbers 
with b < 40 kms -1 and b < 150 kms -1 . For both 
6-samples we also choose three different column den- 
sity ranges: (a) [13.2, 16.5] dex, (b) [13.2, 14.4] dex and 
(c) [14.4, 16.5] dex. The lower limit of sample (a) and 
(b) corresponds to our threshold of completeness. The 
largest observed column density in our sample is ^10 16 ' 5 
cm~ 2 . Hence sample (a) corresponds to the whole sam- 
ple with W > 90 mA, while sample (b) only covers the 
weak Lya systems, and (c) the strong Lya systems. The 
cut at 14.4 was cho sen because the sample analyzed by 
IPenton eT al. (2004) suggested a change of slope near this 
value (see below). Note that for each sample, we only 
consider absorbers with a^/b < 0.4 and cr^/N < 0.4. 

In Fig. [HI we show the column density differential 
distribution of the identified absorbers for the combined 
sample with b < 40 kms -1 and the maximum- likelihood 
fits to the data. The fit to the sample with b < 150 
kms -1 implies a general increase of /3 (see Table [7J for 
the various column density intervals. For sample (a), 
if we vary log A m ; n from 13.2 to 13.4, the results from 
the maximum-likelihood fit essentially do not change. If 
log iVmin < 13.1, (3 decreases rapidly since our sample is 
not complete anymore. If A max is reduced, /3 will in- 
crease as sample (b) shows. For sample (b) with either 
b < 40 kms" or b < 150 kms - if log A m ; n increases by 
up to 0.4 and/or log A max varies by ±J-g dex, the results 
do not statistically change. However, for sample (c), the 
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TABLE 7 

The Power Law Fit to the Distribution 
Function f(N m ) = CmN'f 



[logAr min ,log 


iVmax] 


log Chi 




b < 40 kms -1 


[13.2, 16.5] 
[13.2, 14.4] 
[14.4,16.5] 


1.76 ±0.06 
1.83 ±0.06 
1.52 ±0.10 


12.1 
13.0 
8.4 


b < 150 kms" 1 


[13.2, 16.5] 
[13.2, 14.4] 
[14.4,16.5] 


1.84 ±0.06 
1.92 ± 0.05 
1.61 ±0.11 


13.3 
14.4 
9.9 



Note. — Only data with <j b /b, a N /N < 0.4 are 
included in the various samples (for more details, 
sec $El . 



results are more uncertain. While (3 appears to flatten at 
larger column densities, there are too few absorbers with 
large H I column density to have a full understanding 
of the possible change in the slope. In particular, it is 
not clear from our sample where the break in the slope 
occurs since for log^V m ; n w 14.4 ±0.5 dex, (3 is ~1.5. At 
the 2<7 level, (3 from sample (c) is essentially the same as 
for sample (a) and (b). 

With GHRS and STIS grating moder ate spectral reso- 
lution observations. lPenton et all (|2004l ) found at z < 0.1 
= 1.65 ± 0.07 for 12.3 < logJV(H I) < 14.5 and 
(3 = 1.30 ± 0.30 at log7V(H I) > 14.5 with few data 
points. In their work, 7V(H I) was obtained from the 
equivalent width assuming b — 25 kms~ . Their results 
are consistent w i th th e Key Project dat a presented by 
IWevmann et al.l (|1998[ ) and our results. iDave fc Trippl 
(12001D derived^ = 2.04 ± 0.23 for 13 < log N(R I) < 14. 
They used STIS E140M spectra of two QSOs and were 
able to derive b and N independently using an automated 
Voigt profile fitting software but not allowing for BLAs. 
Within ler, our results are the same. 

At z > 1.5, several studies using high resolution spec- 
tra obtained with the Keck and the VLT have shown 
that (3 ~ 1.5 for H I column density ranging from a 
few times 10 1 2 cm" 2 up to a few times 10 20 cm~ 2 (e.g., 
iTvtler! [19871 : IHu et all [19951 : ILu et alJ fl99fil h although 
there may be so me deviation from a single power law at 
logMH I) > 14 ijPetitiean et alJll993l : IKim et allll997l . 
2002a). There is also some suggestion that for a given 
H I interval, (3 increases as z decreases, but thi s is sta- 
tistically uncertain at z > 2 (IKim et al.ll200"2al ). If we 
compare our results to IKim et all (|2002a[ ) . we find an in- 
crease of (3 as z decreases in the column density interval 
13.2 to 14.4 dex or 13.2 to 16.5 dex, but not in the col- 
umn density interval 14.4 to 16.5 dex, although in this 
range the number of data points is too s mall to draw any 
firm c onclusion. Finally, the analysis of iJanknecht et all 
(2006) at 0.5 < z < 2.0 shows that (3 is intermediate be- 
tween f3 at z < 0.5 and (3 at z > 2. Therefore, it appears 
that the column density distribution steepens with de- 
creasing z in the column density range 13.2 to 16.5 dex. 
A redshift de pendence was f ound in the numerical results 
presented bv lTheuns et alJ (I1998T ). but the observed rate 



6. THE BARYON DENSITY OF THE IGM 

6.1. Narrow Lya Absorption Lines 

To estimate the baryon content of the photoionized 
Lya forest at low z we follow the method presented by 
iSchavdlpOOl . The mean gas density relative to the crit- 
ical density can be obtained from the H I density distri- 
bution function: 

fi(NLA) = >*™*> [ N m ^f(N m )dN m (7) 
p c c J n m 

where mn = 1.6 73 x 10 - 24 g is the atomic mass of hydro- 
gen, /xh = 1.3 corrects for the presence of helium, Hq = 
100 kms^Mpc" 1 , p c = 3H§/(8irG) = 1.06 x 10~ 29 
gcm~ 3 is the current critical density, nu and nm the 
density of the total hydrogen and neutral hydrogen, re- 
spectively, iVm the neutral hydrogen column density, and 
/(Ahi) the differential density distribution function dis- 
cussed in Assuming that the gas is isothermal and 
photoionized, the previous equation can be simplified to, 



fi(NLA) w 2.2 x W- g h 



-l r l/3rp0.59 
I n -LA 



0.16 



1/3 



(8) 



N n 



HI 



of evolution of (3 appears to be smaller than their models 
suggest. 



where h = iV(100kms -1 Mpc -1 ), T 12 is the H I pho- 
toionization rate in units of 10 -12 s , T\ is the IGM 
temperature in units of 10 4 K, and f g is the fraction of 
mass in gas in which s tars and molecules do not con- 
tribute. ISchavd (|20Q1) states that in cold, collapsed 
clumps, f g « 1, but on the scales of interest here f g 
is close to the ratio of the total b aryon density to the 
matter density, which according to ISpergel et al.l (|2003l ) 
is 0.16. We therefore se t f a = 0.16. We use Ti 2 = 0.05 
from lHaardt fc Madaul (fl996h for our average redshift 
z = 0.19. This value of T i 2 is c onsistent with the value 
derived by IDave fc TtTppI |200l[) at z = 0.17. The typ- 
ical temperature of the IGM is not well defined in the 
curren t simulations of the low redshift universe. ISchavd 
(|2001f ) not ed a good a g reeme nt between his estimate of 
T4 and the lDave et alJ (|1999f ) estimate, we therefore use 
Eq.[3]for estimating T4. But as these authors noted there 
is a very large scatter in the T-p/p relation at z ~ 
producing uncertainty in the baryon determination from 
Eq. [5] because of temperature variations. Note for our 
calculation, we set z — z. Because fi(NLA) cx N^~^ 
with (3 > 4/3, the low column density systems dominate 
Q(NLA), therefore we will set the temperature from N m i n 
in Eq. [8] Furthermore since N m { n dominates the solution 
of Eq. [51 using this equation is dependent on the data 
quality of the current sample. We therefore emphasize 
that we only derive the contributions of the denser parts 
of the photoionized Lya forest. 

To estimate the baryon content of the photoionized re- 
gions, we only use H I absorbers with b < 40 kins" 1 . 
Broader systems are here assumed to be principally col- 
lisionally ionized. For the interval of H I column den- 
sity 10 13 2 -10 16 5 cm -2 , using the best fit parameters of 
/(A^hi) listed in Table[7|for this range of column densities 
and with (T 12 , T A , f g , h) = (0.05,1.2,0.16,0.7), we find 
17(NLA)/fi b ~ 0.19 for absorbers with logA^m > 13.2, 
where f^b = 0.044 is the ratio of the total baryon density 
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to the critical density (ISpergel et al.ll2003t iBurles et al.1 
120011 : lO'Meara et~aLl 120011 ). If we assume that the 
power law described in $5] fits the data extending to 
the lower observed column density (12.42 dex), we 
find f2(NLA)/f2b = 0.29 for the column density range 
[12.4,16.5] with (T 12 ,T 4 J g ,h) = (0.05,0.6,0.16,0.7). 
Such an assum ption is not un reali stic since, in high 
redshi ft spectra +Lu et al.1 (|1996l ) and lKirkman fc Tvtlerl 
(|1997h show using their simulation results that weak ab- 
sorbers (down to 12.1-12.5 dex) follow the same column 
density distribution as the stronger absorbers. Assum- 
ing j3 = 1.76 (see Table [7]) and Eq. 02 the dependence 
between f2(NLA)/f2b and A m i n can be approximated by 



ft(NLA)/ft b ~ 0.2 (A min /[10 



13.2 



cm 



2 D 



-0.18 



With the available data, there is no indication of sub- 
component structure in the broad Lya absorption lines. 
However, broadening mechanisms other than thermal 
may also be important. If this is the case, the thermal 
broadening could decrease significantly for the broad Lya 
absorption lines, implying that many of those lines would 
arise in the photoionized IGM, not the WHIM. If the 
complete sample with log Am = [13.2, 16.5] and b < 150 
kms -1 is considered, Q(b < 150kms _1 )/i7b > 0.23 for 
(ri2, T 4 , f g , h) = (0.05, 1.2, 0.16, 0.7). If the sample with 

log A H i = [12.4,16.5] and b < 150 kms" 1 is consid- 
ered, fl{b < 150km s- 1 )/^ = 0.40 for (T 12 ,T 4 J g7 h) = 
(0.05,0.6,0.16,0.7). Therefore, if the BLAs are tracing 
photoionized gas, the estimate of fl would increase by a 
factor ~1.3. These BLAs would not then contribute to 
the baryon budget in the BLAs determined in §6.21 

To estimate a reliable error on fi(NLA) remains a dif- 
ficult task with our current knowledge. The estimate 
of fi(NLA) is model dependent since the representative 
temperature of the IGM is not well known. More numer- 
ical simulations of the low redshift IGM spanning a wider 
range of parameters and using the results from the cur- 
rent sample may tighten L and T. If T changes by 20%, 
J7(NLA)/f2b can change by about 10-15%; if T changes 
by 20%, fi(NLA)/ft b can change by about 10%. The ft 
estimate is also very sensitive to the slope (3: a change 
of (3 by ±0.005 can introduce a change in ft(NLA)/ftb 
by about =f4%. Furthermore the low column density 
absorbers dominate the baryon fraction and the low col- 
umn density cutoff, A m j n , is unknown. Therefore, while 
the exact baryon content is uncertain in the photoion- 
ized IGM, it is clear that the photoionized Lya forest is 
a large reservoir of baryons. The baryon fraction in the 
diffuse photoionized phase traced by (narrow) Lya ab- 
sorbers predicted by cosmological models in the low red- 
shift universe varies by a factor 2 (^2 0-40%) among the 
various simulations (|Dave et al.ll2001f). in general agree- 
ment with our results. We note that lPenton etai] (2004) 
derived ft(Lya)/ftb = 0.29 ± 0.04 for the photoionized 
phase in the column density range [12.5,17.5] dex. 7 In 
view of our discussion on the various uncertainties, we 
believe that their error estimate appears optimistic. 

6.2. Broad Lya Absorption Lines 



7 It is unclear which parameters Pcnton et al. (2004j) used with 
the Schavc (2001) method (if we set the parameters to those given 
byPentonetal.,i.e. (r i2 , T 4 , / 9 , h) = (0.03, 0.5, 0.16, 0.7) we found 
that the values in their Table 4 (Schaye column) are a factor ~3 
too high). 



Baryons also reside in the WHIM, a shock- heated inter- 
galactic gas with temperatures in the range 10 5 to 10 7 K. 
Cosmological hydrodynamical simulations predict that 
the WHIM may contain 30-50% of the baryons at low 
redshift (Cen & Ostriker 1999; Dave et al. 1999). BLAs 
may trace the 10 5 to 10 6 K WHIM if the broadening 
is purely thermal, following T[K] = (6[kms^ 1 ]/0.129) 2 . 
The cosmological mass density of the broad Lya ab- 
sorbers in terms of t oday critical density can be written 
(|Richter et al.ll2004t ISenibach et al.ll2004f ) as, 



ft(BLA) = 



ft(BLA) w 1.667 x 10 



fj,B_mnH T,fn(Tj)N m (i) 
EAA 

SAX 



PcC 



(9) 



where fn(T) being the conversion factor between H I 
and total H and the other symbols have the same mean- 
ing as in Eq. CO For our sample, £ AA = 2.404. In Eq. [9] 
the sum of fn(Ti)N-Ri(i) over index i is a measure of 
the total hydrogen column density in the BLAs. When 
collisional ionization equilibrium (CIE) is assumed, the 
conversio n factor betwe e n H I and total H was approxi- 
mated by iRichter et alj (120041) from the values given in 
ISutherland & Dopital (|1993h for temperatures 10 5 - 10 6 
K: 



log fn(T) « -13.9 + 5.41ogT-0.33(logT) 2 



(10) 



For low-density absorbers, the CIE hypothesis is, 
however, probably a poor approximation because the 
photoionization from the UV ba ckground is important 
(IRichter et al.ll2006bfh Recently. IRichter et all (|2006af) 
us ed the output of the nu merical simulation, carried out 
by iFang fc Brvanl (|2001l ) to investigate the intervening 
O VI absorption in the WHIM, to study the BLAs in the 
low redshift IGM. They found that a significant num- 
ber of BLAs could be photoionized, in particular sys- 
tems that have relatively low 6- values in the range 40-65 
kms -1 . Their simulation that includes collisional ion- 
ization and photoionization suggests that the hydrogen 
ionization fraction can be approximated as: 

log/ H (T)« -0.75+ 1.25 logT. (11) 

According to Fig. 5 in IRichter et alj (|2006aD . the CIE 
hypothesis (Eq. ITU]) provides a firm lower limit to /h (T) . 
Eqs. [10] and [TT] converge as b increases, but never inter- 
sect. 

In Tableland %H1 we saw that the S/N has little 
effect on the estimate for dJV(BLA)/dz for the high- 
est S/N spectra, but can reduce dAf (BLA) / dz for the 
lowest S/N data in our sample. To reduce the uncer- 
tainty from the broadening of the H I absorbers found in 
low S/N spectra, we restricted our sample to absorbers 
with (Tb/b, (Tjv/A < 0.3. To estimate ft(BLA), we con- 
sider the following b- value intervals: [40,150], [40,65], 
and [65,150] kms" . In the [40,65] kms" interval, 
IRichter et al.1 (|2006af ) found that about a 1/3 of the 
BLAs traces cool photoionized gas (T < 2 x 10 4 K), 
i.e. & n t 3> bth for these systems. If all the BLAs are 
assumed to be thermally broadened, we would overesti- 
mate fi(BLA). We therefore randomly pick 2/3 of the 
BLAs in the [40,65] kms -1 &-range to estimate /h- (We 
note that the baryon content of the photoionized BLAs 
with T < 10 4 K is less than 1-2% according to Richter 
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et al. 2006a.) In the [65,150] kms -1 6-range, the num- 
ber of BLAs tracing cool photoionized systems is small 
(less than ~5%) and we assume that all the BLAs in 
this range of 6-values are mostly thermally broadened. 
With the conditions listed above, there are 12 systems 
in the [65,150] kms -1 interval and 31 in the [40,65] 
kms -1 &-range. In Tabled we summarize our estimates 
of fi(BLA)/^ b using Eqs. El QUI and HI] where T was 
de rived assumin g bth ~ 0.96 o bs following the simulation 
of iRichter et al.l (|2006af ). The largest difference in the 
estimates of the baryon budget from CIE and Richtcr 
et al.'s model is for low b- values, where photoionization 
plays a more important role. We note that that if we 
restricted our sample to systems with ab/b, <tn /N < 0.2, 
0(BLA)/51b would decrease by a factor ~1.2, but would 
increase by a factor ^1.2-2.4 if ab/b, <jn /N < 0.4. 

As for the estimate of the baryon budget of the NLAs, 
with our current knowledge it appears a very difficult 
task to estimate an error on the baryon budge t in BLAs. 
It is model dependent. IRichter et alj (|2006al ) provided 
estimates of the amount of the BLAs that probe cool 
photoionized gas and the amount of thermal broaden- 
ing in BLAs. These estimates may change with refined 
cosmological models. If we assume that the broaden- 
ing for all the BLAs is purely thermal (b t h — b b s ), 
the estimates of tt(BLA)/tt b in the [40,150] kms -1 b- 
range would increase by a factor 1.6, i.e. in the CIE 
model, f2(BLA)/J7b — 0.13, and i n the collisional ioniza - 
tion plus photoionization model of lRichter et al.l ((2006a) , 
fi(BLA)/^b — 0.32. The broadening could be due to 
several unresolved components in some cases or the non- 
thermal broadening could be important. This is partly 
accounted for since we remove a 1/3 of the BLAs in 
the [40,65] kms -1 &-range. Since we discarded sys- 
tems with large errors in b and N, for the current data 
set, a single Gaussian fit provides a good representa- 
tion of the observed profiles. However, as we discussed 
in Q3.2\ the formal errors given by the fit may not be 
sufficient especially f or the weaker and broader systems. 
IRichter et al.l (|2006bl ) and the other papers mentioned in 
Table [T] present several examples of BLAs that are well 
fitted with a single Gaussian with sometimes the presence 
of several H I Lyman series transitions that further con- 
strain the fit parameters (see also the discussion in §4.2j) . 
To further explore the quality of the profile fitting, the 
Cosmic Origin Spectrograph (COS) could really improve 
the situation by increasing the number of sightlines and 
increasing the S /N in the spectra. Additional theoretical 
simulations of the low-z IGM with higher resolution and 
refined physics combined with the analysis of simulated 
spectra should improve our understanding of the BLAs. 

In summary, the CIE estimate provides a firm and 
conservative limit for the baryon budget in collision- 
ally ionized BLAs, f2(BLA)/J7b > 0.1 for systems with 
40 < b < 150 kms -1 , logiVm > 13.2, <r b /b,cr N /N < 0.3. 
The BLA simulations show that the ionization fraction 
in BLAs is governed by collisional ionization and pho- 
toionization. Therefore, il(BLA)/J7b > 0.2 (with the 
same constraints that those listed above) should reflect 
better the baryon budget of the BLAs in the low-z IGM. 
This is likely a lower estimate as well since the existing 
observations do not have the S/N to detect the broader 
systems (where most of the baryons are believed to exist) 



TABLE 8 

Estimates of the Baryon Density in 
BLAs 



b [kms- 1 ] 


n CIE (BLA) a 


n c+P (BLA) 6 


n b 


n b 


[40, 65] 


0.03 


0.11 


[65, 150] 


0.05 


0.09 


[40, 150] 


0.08 


0.20 



Note. — Only data with a b /b, a N /N < 0.3 
and log -/Vhi > 13.2 are included in the vari- 
ous samples. In the [40,65] kms -1 interval, 
only 2/3 of the BLAs are considered. We de- 
rive the temperature of the absorbers assum- 
ing 6th — 0.96 o bs- The baryon density in the 
interval [40,150] kms" 1 is the sum ofH(BLA) 
in the two other intervals, a: Obtained from 
Eqs. [9] and |10| baryonic density for the BLAs 
assuming that they are collisionally ionized 
and CIE applies, b: Obtained from Eqs.[9]and 

1 1 X | baryonic density for BLAs assuming that 
they are collisionally ionized and photoionized 
using the res ults from the hydro dynamical 
simulations of lRichter et all j2006aD (see $57Z\ 
for more details). £2 b = 0.044 JSpergel et al.l 

120031) . 

and the numerous weaker BLAs (see Fig. EJ) ? an d are not 
complete at b > 80 kms - . 

6.3. Baryon Budget of the Low Redshift Universe 

At redshift 2-4, observations of the Lya forest show 
that about 80-90% of the total baryo n budget is found in 
the cool photoionized phase of IGM (TRauch et al.|[T997l ; 
IWeinberg et al.lll997t Ikim et al.ll200lD . In the low red- 
shift universe, " iFukugita et all (|1998f) show that 9-10 % 
of the total baryons are found in galaxies (stars, neutral 
gas, molecular gas) and ~6 % of the baryons are associ- 
ated with the hot plasma in clusters of galaxies. 

Our analysis of the low redshift Lya forest shows the 
presence of NLAs and BLAs. The NLAs trace the cool 
photoionized IGM, while the BLAs are likely to probe the 
cool photoionized IGM and the hot highly-ionized gas in 
the WHIM. There are still many uncertainties that are 
not controlled in the determination of the baryon content 
in the low redshift Lya forest. But at least 20% of the 
baryons are found in the denser parts of the photoion- 
ized phase (NLAs, b < 40 kms -1 , log Nm > 13.2, and 
ab/b, vn/N < 0.4) and at least ~10% are found in the 
WHIM (BLAs, 40 < b < 150 kms -1 , \ogN m > 13.2, 
ab/b, aN/N < 0.3), corresponding to at least 30% of the 
total baryon budget. We believe that these limits are 
very conservative since they do not take into account 
the low column density photoionized absorbers and the 
broadest BLAs. Furthermore recent simulations of the 
BLAs in low redshift universe show that the BLAs may 
trace both photoioniz ed and collisionally ion ized gas. Us- 
ing the simulations of lRichter et al. (2006aj), we find that 
the BLAs include at least 20% of the baryons. We have 
also seen that fi in the cool photoionized IGM is domi- 
nated by the low column density H I absorbers (i.e. ab- 
sorbing gas with overdensity p/p < 2). If the weak sys- 
tems follow the same differential column density distri- 
bution function as the stronger absorbers, the cool pho- 
toionized gas traced by the NLAs with logA^i > 12.4 
can contain ~30% of the total baryons. Combining these 
NLA and BLA budgets, the low redshift IGM contains 
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at least ^50% of the baryons. This is much larger than 
the known amount of baryons in galaxies and clusters 
of galaxies in the low redshift universe. High S/N data 
will allow us to search for the weakest Lya absorbers 
(< 10 125 cm~ 2 ) and estimate their column density dis- 
tribution, search for the very broad absorbers and the 
weak BLAs, and provide a better understanding of their 
physical nature. It is, however, already apparent from 
this study that the IGM traced by the narrow and broad 
Lya absorption lines in the low redshift universe is a ma- 
jor reservoir of baryons. 

7. SUMMARY 

We analyze the physical properties of the low-z IGM 
traced by H I absorbers from a sample of 7 QSOs ob- 
served with the E140M mode of STIS (R ~ 44000) and 
with FUSE (R ~ 15000). These seven lines of sight were 
fully analyzed and presented in recent and future papers 
(see references in Table Q] and the Appendix of this pa- 
per). Our sample has a total unblocked redshift path 
of 2.064 corresponding to a total absorption distance of 
2.404 and is complete for log Ahi > 13.2 provided b < 80 
kms -1 . With high spectral resolution the column den- 
sity (JVhi) and the Doppler parameter (b) were indepen- 
dently estimated, providing the opportunity to directly 
study the distribution and evolution of b. The spectral 
resolution of the STIS E140M H I observations in our 
sample is similar to that of the high redshift (z > 1.5) ob- 
servations obtained with Keck/HIRES and VLT/UVES. 
This similarity allows a relatively simple study of the 
evolution of the 6-parametcr of the neutral hydrogen gas 
in the intergalactic medium. Because the nominal lower 
temperature of the WHIM is T ~ 10 5 K, we consider 
two different populations of H I absorbers throughout 
this work: narrow Lya absorbers (NLAs) have b < 40 
kms~ and the broad Lya absorbers (BLAs) have b > 40 
kms~ . Our main findings applicable at redshifts z < 0.4 
are summarized as follows: 

1. The H I Doppler width distribution has a high ve- 
locity tail that develops at b ~ 40-60 kms -1 . The b-Nm 
distribution reveals for the NLAs an increase of b with in- 
creasing Ahi with a very large scatter. We find that most 
of the BLAs are found in the H I column density range 
[13.2,14.0] dex and the broader systems are found for 
13.1 < log Am < 13.5. Recent cosmological simulations 
of BLAs show a good agreement with the observations, 
but clearly indicate that broader absorbers and many 
weak BLAs (log Ahi < 13.1) remain to be discovered. 

2. We find dN (BL A) / 'dz = 30 ± 4 for absorbers with 
40 < b< 150 kms" 1 , logA ffl > 13.2, and a b /b,a N /N < 
0.4. The narrow H I absorbers are more frequent, with 
dA/"(NLA)/dz = 66 ±6 for absorbers with b < 40 kms -1 , 
log Am > 13.2, and a^/b^N /N < 0.4. Very narrow 
absorbers (b < 15 kms" 1 ) with log Am > 13.2 are scarce. 
The number of weak NLAs and BLAs (log A H i < 14) is 
far larger than the number of strong NLAs and BLAs, 
respectively. 

3. For the narrow absorbers with 13.2 < log Am < 
16.5 and Uhjb.u^ /N < 0.4, we find that the column 

density distribution, /(Am) oc N^f , can be fitted with 
(3 = 1.76 ± 0.06. For the entire sample (i.e. including 
NLAs and BLAs), the slope changes to 1.84 for the same 
column density range. We confirm an increase of p with 



decreasing z when our results are compared to higher 
redshift analyses. There is some weak evidence for a 
break at log Am ~ 14.4, but the location of this break is 
uncertain because of the small number statistics for the 
higher column density lines. 

4. We argue that various samples probing different 
redshift ranges can be directly compared when the con- 
ditions < b < 100 kms" 1 , 13.2 < logA ffl < 14.0, and 
crf,/6,CTjv/A < 0.4 are set. The distribution of b for the 
broad absorbers has a distinctly more prominent high 
velocity tail at low z than at high z (1.5 < z < 3.6). 
The median and mean b- values are systematically larger 
by -15-30% at z < 0.5 than at 1.5 < z < 3.6. The 
ratio of the number density of BLAs to NLAs at low z is 
larger than at high z by a factor —3. This suggests that 
a larger fraction of the low-z universe is hotter than at 
1-5 < z < 3.6 and/or the low- z universe is more kine- 
matically disturbed than the high-z universe. 

5. The NLAs trace the cool photoionized IGM 
at T < 10 4 K, and we find that Q{NLA)/n h ~ 
0.2(AW[10 13 - 2 cm- 2 ])- - 18 (when f(N m ) cx A^' 76 
and assuming Am oc T~ 0,42 ), where A m ; n is the low- 
est H I column density in the sample considered to esti- 
mate fi(NLA). The contribution to the baryon budget 
of the denser parts of the photoionized Lya forest with 
log N m > 13.2 is ~20%. If the weakest H I absorber with 
logAm = 12.4 follows the same /(Am), the cool pho- 
toionized gas with log Am > 12.4 contains about 30% of 
the total baryons. 

6. The BLAs trace mostly the highly-ionized gas in 
the WHIM at T ~ 10 5 -10 6 K if the width of the BLAs 
is dominated by thermal motions. The assumption of 
collisional ionization equilibrium provides a firm lower 
limit to the amount of the baryons in the WHIM traced 
by the BLAs (40 < b < 150 kms" 1 , logA ffl > 13.2, 
a b /b,a N /N < 0.3): tt(BLA)/n h > 10%. If a hydro- 
dynamical simulation including the effects of collisional 
ionization and photoionization is used to estimate the 
ionization correction for the BLAs, we find that at least 
20% of the baryons are in the BLAs. 

7. Our most conservative estimate of the baryon 
budget in absorbers with logAni > 13.2 shows that 
the Lya forest has at least about 30% of the baryons 
(NLAs+BLAs), far larger than the baryon budget in 
galaxies. We suggest that the low redshift IGM with 
T < 10 6 K traced by the NLAs with logA H i > 12.4 
and the BLAs with log Ahi > 13.2 could contain at least 
50% of the baryons (this is a lower limit since the broader 
BLAs where most the baryons reside remain to be dis- 
covered). The estimates of the amount of baryons still 
rely, however, on critical assumptions (e.g. pure thermal 
broadening versus other broadening mechanisms, tem- 
perature of the IGM, behavior of the weak narrow sys- 
tems) that produce uncertainties that are currently not 
well controlled. Observationally, if COS is deployed, our 
understanding of the BLAs, in particular, will improve 
given the high S/N spectra that COS should obtain (but 
with approximately 2 times lower resolution than for the 
STIS observations presented here). The present paper 
and future observations of the low-z IGM should moti- 
vate more precise cosmological simulations of the low- 
and high-z IGM that are needed for our understanding 
of the intrinsic properties of the BLAs and NLAs. 
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APPENDIX 

We have revisited the H I measurements presented by IWilliger et all (|2006( ) toward PKS 0405-123. We were first 
motiv ated to produce a new analysis of the PKS 0405-123 FUSE and STIS/E140M spectra because IWilliger et all 
( 200|) noted that several of their BLAs may be noise features. For our work, we needed to clearly differentiate a real 
detection from a noise feature. While reanalyzing the BLAs, we noted that some narrow H I lines were either not 
present in the spectra at the 3<r level (although Williger et al. claimed to report only Aa features) or were misidcntificd. 
We have not identified the ultimate origin of the problem, but we suspect that the difference with Williger et al. can 
ultimately be traced to a difference in the manner in which the continuum placement was done. Williger et al. fitted a 
single global continuum, while we determine continua separately placed within ±1000-3000 kms -1 of each absorption 
line. We therefore produce here a new line list and new measurements for H I toward PKS 0405-0123 in Table Al. 
This reanalysis also provides an overall coherent data sample, since the other sightlines were analyzed following the 
same methodology. In Table Al, we report our measurements for H I toward PKS 0405-123. These H I absorbers 



Lya Forest at z < 0.4 



25 



are detected at least at the 3tr level. We list in this table the observed wavelengths and the rest-frame equivalent 
widths of the Lya transition. The listed column dens ities and Doppler p arameters were measured in the rest-frame 
by N. Lehner (NL) following the method described in lLehner et al.l (|2006f ). The weak absorbe r s (log iVHi < 14) were 
also independently measured by P. Richter (PR) using the method described in iRichter et al.l (|2004T ). In particular, 
different profile fitting software wa s used by NL and PR . In Fig.ll5( we compare the column densities and b- values (for 
the weak absorbers) presented by IWilliger et al.l (|2006f ) to those listed in Table Al (left panels). On the right-hand 
side, we present the measurements obtained by PR against those listed in Table Al. For both b and Nm, there is 
an excellent agreement between PR and NL measurements at the la level. In contrast, Williger et al. (GW) found 
generally larger column densities and Doppler parameters. The main difference between our analysis (NL and PR) and 
GW's analysis is that we systematically and independently chose the continua within ±1000-3000 kms - of the H I 
absorption line, while GW fitted the continuum using a semi-automatic procedure (AUTOVP) combined with custom 
fitting in the vicinity of the QSO emission lines. Since these weak absorbers are generally detected solely in Lya, the 
continuum placement is critical. Our method for the continuum placement was followed for all the other sightlincs 
presented in this work. For the absorbers where Lya and Ly/3 are detected, our measurements and those made by 
GW agree within la. 

Absorbers with logAni > 14 are detected in several Lyman series lines. For the absorbers that were fitted with 
a single component (z = 0.03000,0.35099,0.40571,0.40886), GW and our measurements are in agreement within la. 
GW fitted the systems listed in Table Al at z = (0.09657,0.09659) and (0.18262,0.18287) with a single component. 
We found that a 2-component fit provides a statistically better fit. The total column densities for those absorbers are 
in agreement within la with those measured by GW. For the partial Lyman limit system at z = 0.167, we found that 
5 components provide a better fit to this complicated absorber (see Fig. [T^l and note (6) in Table Al). We note that 
our fit is not yet as good as it should be despite our best effort (see Fig. [16]). In particular, the broad component 
at z = 0.16678 may actually be narrower since metals- ions such as C III are detected at this redshift (this directly 
affects the measurements of the absorber at z — 0.16661 where C III is also detected). Finally, we fitted the system 
at z — 0.36808 with 4 components. This is not statistically different from the 2 component fit presented by GW, but 
the 4-component fit appears to be a better model when the kinematics of the metal-ions are taken into account (see 
note (12) in Table Al for more details). 

In Table A2, we report the list of features that were identified as Lya and claimed to be detected at least at the 
4a level by GW. These features are actually either non-detections at the 3tr level (according to our measurements) or 
were misidentified. The 181 Lya absorbers toward PKS 0405-123 reported by GW now reduce to 74 absorbers (and 
in a few cases, those are multi-component absorbers). 



Lehner et al. 



TABLE 9 

IGM H I REVISED MEASUREMENTS TOWARD PKS 0405-123° 
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TABLE 9 — Continued 



X 


z 




log JV H i 


6hi Note 


(A) 




(ml) 


(dex) 


(kms" 1 ) 



NOTE. — a: The wavelength, A, is the observed redshifted wavelength of Lya. The equivalent width, WL ya , is the rest-frame equivalent width 
of Lya. A colon means that the result is uncertain. 

(1) The Lya for the systems at z — 0.05925 and 0.05896 were fitted simultaneously. 

(2) We fitted simultaneously detected H I AA1215, 1025, 972, 937 (H I A949 is blended). 

(3) We fitted simultaneously detected H I AA1215, 1025, 972, 937 (H I AA972, 949 are blended). We found that a 2-componcnt fit improves the 
modeling of Lya, yielding a nar row and broad compon ent. If the BLA is mostly thermally broadened, it yields a temperature of T ~ 3 X 10 J K. 
O VI is found in this absorber (Prochaska ct al. 2004). The temperature found from the H I broadening is consistent with O VI being mainly 
collisionally ionized. 

(4) This profile appears asymmetric, but within the S/N a single component fit is suitable. This system is, however, not accounted in the BLAs 
in this paper because of the asymmetry. 

(5) The Lya for the systems at z — 0.16121 and 0.16150 were fitted simultaneously. 

(6) The systems at z = 0.16661, 0.16678, 0.16694, 0.16712, and 0.16714 were fitted simultaneously. The equivalent width re ported at z = 0.16714 is 
the total equivalent width for these absorbers. The systems at z — 0.1671 are detected from Lya down to th e Lyman Limit. [Prochaska ct al. (2004) 
derived the total column density of the system at z ~ 0.1671 using the flux decrement at the Lyman limit. Willigcr ct al. (2006) fixed the column 
density found by Prochaska ct al. (2004) for the systems at z — 0.16714 and used H I AA1215, 1025 to derive b. They fitted a very broad absorber 
at z — 0.16692 and a very narrow absorber at z — 0.16628. This fit appears satisfactory for H I AA1215, 1025, but not for the other transitions (see 
Fig. 1 16 1 ■ In particular, the core of the higher H I Lyman scries is not reproduced. To reproduce the core of the line, we find that a 2-componcnt fit 
is necessary: a very narrow component with b — 12 kms -1 and a broader (but still relatively narrow) component with 6—30 kms" 1 . Furthermore 
the C II and C III kinematics (see Figs. 9 and 10 in Prochaska et al. 2004) indicate 4 distinct components at —136, —93, —40, kms" 1 . Wc, 
therefore, fitted 5 components set initially at —136, —93, —40, —5, 5 kms -1 . The velocity of the component at —40 kms -1 was not allowed to 
vary; all the other parameters for the various components were allowed to vary. For the fit, we used simultaneously detected H I AA1215, 1025, 
937, 930, 926, 923, 919, 918, 917 lines. H I AA972 and 949 arc contaminated, while H I A920 suffers from fixed-pattern noise. Wc did not use 
higher Lyman seri es transitions because t he continuum placement becomes uncertain due to the flux decrement. The resulting fit is shown in Fig. 

mWc note that Prochaska ct al. (2004) used a curve-of-growth analysis with two components separated by 40 kms based on the C II and 
I profiles and found 16.35 and 15.5 dex, which is ro ughly consistent with our results for the systems z — 0.16694,0.16712,0.16714. While our 
fit appears more satisfactory than the one presented by Williecr ct al. ( 2006) who did not take into account the velocity structure of the metal-ion 
profiles, the solution may actually be more complicated. The errors do not reflect that there may be more than 5 components. In particular, the 
broad component at z — 0.16678 may actually be narrower since metal-ions such as C III are detected at this redshift; this directly affects the 
measurements of the absorber at z — 0.16661 where C III is also detected. Therefore we mark the measurements at z — 0.16661 and 0.16678 as 
uncertain. 

(7) The syste ms at z = 0.18261 and 0.18287 were fitted simultaneously. We used detected H I AA1215, 972, 949 (H I AA1025, 937 are contaminated). 
Wc note that Willigcr ct al. (2006) fitted a single component at z — 0.18271: a 2-componcnt fit improves the fit. The equivalent width is the total 
equivalent width including both systems. 

(8) The systems at z — 0.24513 and 0.24554 were fitted simultaneously. Wc used detected H I AA1215,1025 in the profile fitting. 

(9) Wc fitted simultaneously detected H I AA1215, 1025. Wc note that H I A972 is also detected and possibly contaminated by a weak Lya system. 

(10) We fitted simultaneously detected H I AA1215, 1025. 

(11) The Lya for the systems at z — 0.34186 and 0.34234 were fitted simultaneously. 

(12) The equivalent width is the total equivalent width including the systems at z — 0.36079, 0.36080, 0.36115, 0.36150. The system at z — 0.36080 
is detected in H I AA1215, 1025, 972, 949, 937, 930, 923 (H I A926 is contaminated). We undertook to fit these lines with 2, 3, and 4 components. 
Statistically, these fits are comparably good because Lya is so noisy. The 3- and 4-componcnt fits appear, however, to be better models than the 
2-componcnt fit, especially when one takes into account that O VI AA1031, 1037 appear to be present in the weakest component. The 3-componcnt 
fit docs not fit the blue edge of the Lya profile as well as the 4-componcnt fit, but other wise the 3- or 4-comp onent fits are quite similar. We note 
that our 2-component fit yields similar results to those found by the 2-component fit of[Williscr ct al. ( 2 00611 . 

(13) Partially blended with interstellar C I* A1657.379. 

(14) We fitted simultaneously detected H I AA1215, 1025, 972, 949, 937, 930, 926, 923. 

(15) We fitted simultaneously detected H I AA1215, 1025, 972, 949. 

TABLE 10 

NON- DETECTION AN D CONFUSION WITH OTHER ISM/IGM LINES LISTED AS 
Lycv DETECTION IN IWlLLIGER ET AlI [|2006D TOWARD PKS 0405-123° 
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giVfil 


Note 
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1218.8185 
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1220.3381 
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A 2 Wm logJVm Note 

(A) (mA) (dex) 
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.81 


1465. 


.3808 





.20541 


< 


36 


,6 


< 


12 


.83 


1467. 


.3259 





.20701 


< 


32 


.6 


< 


12 


.77 


1467. 


.8486 


o 


.20744 


< 


33 


,4 


< 


12 


.79 


1474 


5713 





21297 


< 


42. 


.4 


< 


12 


.89 


1474. 


9602 





21329 


< 


42 


.3 


< 


12 


.89 


1476. 


.2732 





21437 


< 


42 


.6 


< 


12 


.89 


1477. 


0876 


o 


.21504 


< 


34 


.6 


< 


12 


80 


1478 


.8383 


o 


.21648 


< 


37 


5 


< 


12 


.84 


1479. 


.5311 





.21705 


< 


36 


,2 


< 


12 


.82 


1480 


0175 





21745 


< 


35. 


.7 


< 


12 


.82 


1480 


6374 


o 


.21796 


< 


34. 


4 


< 


12 


80 


1481. 


2574 


o 


.21847 


< 


35 


6 




12 


82 


1482. 


.9472 


o 


21986 


< 


35 


.8 


< 


12 


.82 


1484 


.4182 





.22107 


< 


38. 


,7 


< 


12 


.85 


1487. 


.4573 





.22357 


< 


39 


.6 


< 


12 


.86 


1488. 


.6122 


o 


.22452 


< 


39 


,4 


< 


12 


.86 


1489. 


.4268 





.22519 


< 


39 


6 


< 


12 


.86 


1490. 


0102 





22567 


< 


39 


3 


< 


12 


.86 


1493. 


5965 





22862 


< 


38 


,7 


< 


12 


.85 


1494. 


.5447 





22940 


< 


38. 


,7 


< 


12 


.85 


1495 


1282 


o 


.22988 


< 


39 


1 


< 


12 


86 


1499. 


8450 





.23376 


< 


39 


.6 


< 


12 


.86 


1501 


.5835 





.23519 


< 


39 


2 


< 


12 


.86 


1506 


.5190 





23925 


< 


40. 


1 


< 


12 


.87 


1507. 


.1147 


o 


.23974 


< 


40. 


.9 


< 


12 


.88 


1510. 


.3605 


o 


24241 


< 


39 


.6 


< 


12 


86 


1510 


7740 


o 


.24275 


< 


39. 


4 


< 


12 


86 


1517. 


.8612 





.24858 


< 


52 


.5 


< 


12 


.99 


1538. 


.0050 





.26515 


< 


43 


,7 


< 


12 


.91 


1541 


2264 





.26780 


< 


45. 


,6 


< 


12 


.92 


1543 


0499 


o 


.26930 














1544. 


.3142 





27034 


< 


45 


.1 


< 


12 


.92 


1544 


.9585 





.27087 


< 


45. 


1 


< 


12 


92 


1545 


5785 





.27138 


< 


42. 


,7 


< 


12 


90 


1561. 


.5525 





.28452 


< 


47 





< 


12 


.94 


1568. 


.7249 





.29042 


< 


49 


,7 


< 


12 


.96 


1569. 


.4422 





29101 


< 


49. 


.1 


< 


12 


.96 


1573. 


3566 





.29423 


< 


42 


.5 


< 


12 


.89 


1580. 


3102 





.29995 


< 


57 


.4 


< 


13 


.02 


1586 


.4007 





.30496 


< 


57. 


4 


< 


13. 


02 


1588. 


.5768 





.30675 


< 


47. 


,5 


< 


12. 


,94 


1589. 


.1116 





30719 


< 


45 


.5 


< 


12 


,92 


1592 


.3453 





30985 


< 


48 


9 


< 


12. 


95 


1594 


.6187 





31172 


< 


48. 


9 


< 


12. 


,95 


1595. 


.1779 





31218 


< 


50 


.1 


< 


12 


,96 


1596. 


.1139 





31295 














1596. 


.5880 





.31334 














1597. 


.1716 





.31382 














1606. 


0945 





.32116 


< 


48 


.5 


< 


12 


,95 


1614. 


.9324 





32843 


< 


62 





< 


13. 


06 


1615 


.6011 





.32898 


< 


59. 


9 


< 


13. 


,04 



3 



Lya Forest at z < 0.4 29 



TABLE 10 — Continued 



A 




z 


Whi 


log iVm Note 


rAi 






(mA) 


(dex) 


1623.9163 





33582 


< 44.5 


< 12.91 


1624.9982 





33671 


< 45.7 


< 12.92 


1625.6182 





33722 


< 44.9 


< 12.92 


1628.5358 





33962 


< 46.2 


< 12.93 


1636.8388 





34645 


< 52.5 


< 12.99 


1639.2459 





34843 


< 73.5 


< 13.13 


1668.3856 





37240 


< 61.3 


< 13.05 


1673.1509 





37632 


< 83.0 


< 13.18 


1678.6093 





38081 


< 51.8 


< 12.99 


1684.1528 





38537 


< 49.3 


< 12.96 


1685.4779 





38646 


< 50.4 


< 12.97 


1696.8565 





39582 


< 64.2 


< 13.07 


1700.7588 





39903 


< 62.2 


< 13.06 


1702.0596 





40010 


< 62.0 


< 13.06 



Note. — a: The upper limits are 3cr and were estimated over a velocity range [—40,40] km s 

(1) Misidentification: interstellar C I* A1329. 

(2) The spectrum is erratic near this wavelength. 

(3) Misidentification: O VI A1031 at z = 0.4951. 

(4) There is undulation in the spectrum, making any identification with a BLA very uncertain. 
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Fig. 15. — Comparison of H I column densities and Dopplcr parameters for the weak absorbers (log./VHi < 14). The left-hand side 
compares the fit results from Willigcr et al. (2006, GW) and N. Lchncr (NL). The right-hand side compares the independent measurements 
produced by P. Richter (PR) and N. Lehner (NL). Identical results would lie along the solid line. 
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Fig. 16. — Profile fittin g to the system at z = 0.167. The solid thick line shows the fit adopted in Table Al. The dotted thick line 
shows the fit adopted by Willigcr ct al. (2006). The light parts of the spectra denote blends with other features. The tick marks in the 
H I A1215 panel are the centroids of the 5 components found in our profile fitting. From left to right, the tick marks correspond to 
z = 0.16661, 0.16678, 0.16694, 0.16712, 0.16714. Discrepancies of the fits to the data are discussed in the text and footnote 6 of Tablc[9] 



